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ABSTRACT 


Analytical  results  for  estimating  natural  frequencies  and  stress  response  for  conventional 
stiffened  flat  panels,  two-sided  box  structure,  wedge  structure,  and  unstiffened  cylin¬ 
drical  panels  are  presented.  Acoustic  fatigue  tests  of  12  stiffened  flat  panel  and  six  box 
structure  designs,  two  specimens  for  each  design,  were  used  to  develop  design  equations 
for  estimating  acoustic  fatigue  resistance  of  skin,  stringer,  and  rib  structure.  Computer 
programs  are  presented  for  estimating  frequencies  of  stiffened  flat  panels,  box  structure, 
and  wedge  structure.  Nomographs  are  presented  for  skin  design  and  cur'-ed  panel  stress 
response . 
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I  -  INTRODUCTION 


Previous  investigations  directed  toward  eliminating  or  reducing  acoustic  fatigue  effects 
from  intense  noise  on  flight  vehicle  structure  have  established  design  criteria  and  nomo¬ 
graphs.  The  original  development  program  was  conducted  by  McGowan  (ASD-TDR- 
63-820,  1963).  Nomographs  were  empirically  derived  from  acoustic  fatigue  tests  of 
certain  popular  structural  designs.  The  fatigue  tests  were  conducted  generally  using 
discrete  frequency  excitation,  and  the  resulting  fatigue  data  were  converted  to  'random 
fatigue  data1  through  the  use  of  Mile's  single  degree-of-freedom  theory  and  the  Miner- 
Palmgren  cumulative  damage  rule.  These  nomographs  have  been  helpful  in  designing 
acoustic -fatigue-resistant  structure  for  present-day  and  near-future  aircraft. 

For  high-performance  aircraft,  lightweight,  acoustic -fatigue-resistant  structure  is 
essential.  To  meet  these  critical  demands  for  weight  reduction  and  increased  resistance  to 
acoustic  fatigue,  McGowan's  work  was  refined  by  Ballentine  (AFFDL-TR-67-156,  1968) 
with  the  objective  of  removing  conservatism  and  extending  the  range  of  application  of  the 
existing  design  methods.  Ballentine  conducted  additional  acoustic  fatigue  tests  using 
random  noise  excitation.  The  structural  configurations  considered  by  Ballentine  consisted 
of  flat  stiffened  panels  and  flat  honeycomb  sandwich  panels. 

McGowan's  original  work  concerning  acoustic  fatigue  resistance  of  internal  airframe 
structure  (stringers,  ribs,  and  frames)  has  proved,  in  practice,  to  be  conservative  from  a 
weight  standpoint.  Ballentine  based  the  design  of  his  stiffened  panel  specimens  on 
McGowan's  design  techniques.  Of  the  sixty  stiffened  panel  specimens  tested  by 
Ballentine,  no  failure  of  the  internal  structure  was  observed.  Hence,  an  exploratory 
program  to  define  more  accurately  the  acoustic  fatigue  resistance  of  internal  structure  was 
required  in  order  to  remove  the  weight  penalty  incurred  by  using  available  design 
techniques . 

Bailentine's  work  was  the  first  step  in  the  acoustic  fatigue  refinement  program.  The 
research  study  described  here  is  a  continuation  of  this  refinement  program.  A  description 
of  the  analytical  investigation  is  presented  in  Section  II;  the  experimental  program  is 
described  in  Section  III;  Section  IV  presents  the  correlation  of  test  data  and  analytical 
results;  and  the  development  of  design  equations  and  nomographs  is  presented  in  Section  V. 
Appendices  are  included  that  describe  the  data  collection  and  reduction  system,  detailed 
specimen  mode  shapes,  and  computer  programs  for  calculating  frequencies  and  mode  shapes 
of  stiffened  flat  panel  structure,  two-sided  structure,  and  wedge  structure. 
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II  -ANALYTICAL 


A.  Introduction 

The  main  emphasis  of  this  program  was  to  provide  experimental  data  for  the  acoustic  fatigue 
resistance  of  internal  airframe  structure.  The  analytical  effort  was  conducted  to  suggest 
the  significant  parameters  and  mathematical  expressions  governing  orthogonally  stiffened 
flat  panels,  box  structure  (typical  of  flaps),  wedge  structure  (typical  of  control  surfaces), 
and  simple  cylindrical  panels.  The  structural  configurations  considered  analytically  were 
selected  based  upon  the  experimental  specimen  designs.  In  order  to  develop  simple 
expressions  suitable  for  design  application,  several  approximations  have  been  made  in 
formulating  the  analytical  models.  These  approximations  allow  the  development  of  design 
equations  for  estimating  the  significant  parameters  in  a  form  suitable  for  design  studies. 


B.  General  Theory 

1 .  Response  of  Linear  Sing le-Degree-of-Freedom  Oscillator 
to  Random  Acoustic  Excitation 


In  Reference  1,  Clarkson  presents  a  summary  of  the  general  theory  of  response  of  a  linear 
oscillator  to  random  acoustic  excitation  .  The  bas  e  assumptions  in  the  development  are 
that  only  one  mode  is  predominant  in  the  frequency  range  of  interest,  that  the  excitation 
spectrum  changes  slowly  with  frequency  near  the  predominant  mode,  and  that  the  excita¬ 
tion  pressure  is  completely  in  phase  over  the  portion  of  structure  of  interest.  For  t_hese 
assumptions,  the  mean  square  displacement  response  for  the  r^  mode  at  the  point  x  on  the 
structure  is  given  by  the  equation 

2 

G  (u)  (1) 

In  terms  of  the  static  displacement  W  (x)  due  to  a  uniform  pressure  of  unit  magnitude,  the 
mean  square  displacement  can  be  wriften  as 
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where 
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Fora  linear  single-degree-of-freedom  system.  Miles  has  shown  in  Reference  2  that  the 
mean  square  stress  at  the  point  of  interest  can  be  expressed  as 


"J  ■'  '””r 
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(3) 
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where  the  viscous  damping  ratio,  Q  ,  and  the  response  frequency,  f  (Hz),  for  the  r  mode 
have  been  introduced .  The  term  a  is  the  dimensionless  stress  at  the  point  of  interest  due  to 
a  uniform  static  pressure  of  unit  magnitude.  The  emphasis  of  the  analytical  program  is  to 
develop  expressions  for  the  natural  frequency,  f  ,  of  the  r™  structural  mode  and  the  dimen¬ 
sionless  stress,  a  ,  as  indicated  in  Equation  3  . 


2,  Frequency  Analysis 

Single  mode  response  of  a  structure  can  be  estimated  using  the  Rayleigh  method  provided 
that  reasonably  accurate  mode  shapes  are  used  in  the  analysis.  For  the  structural  configu¬ 
rations  considered,  the  Rayleigh  energy  method  has  been  used  exclusively.  The  main 
emphasis  was  to  develop  suitable  mode  functions  for  the  structure.  Letting  Wr  denote  the 
amplitude  of  the  r™  mode  of  the  structure,  the  kinetic  energy  and  the  potential  energy  can 
be  expressed  as 

T=Jm  w2  U  =^K  W2  (4) 

2  r  r  2  r  r 


where  M  is  the  moda  I  mass . 
r 

K  is  the  modal  stiffness, 
r 

For  simple  harmonic  motion,  the  kinetic  energy  can  be  expressed  in  terms  of  the  radian 
frequency,  ti^,  as 


T  = 


1  2  2 
%  M  VT 
2  r  r  r 


(5) 


Equating  the  maximum  values  of  the  kinetic  and  potential  energy  one  obtains  for  simple 
harmonic  motion  the  expression 


u)r=V(K  (6) 

The  definitions  in  Equations  4  and  6  will  be  used  to  estimate  the  natural  frequency, 
f  =i«  /2tt,  for  the  structure  and  mode  under  consideration  .  Expressions  for  the  dimen¬ 
sionless  stress,  Oq,  will  be  developed  subsequently  for  each  structure  considered. 

C.  Flat  Stiffened  Panels 

The  structural  configuration  under  consideration  is  illustrated  in  Figure  1  .  The  nine-bay 
panel  is  typical  of  configurations  adopted  for  acoustic  fatigue  testing.  The  portion  of  the 
structure  of  interest  is  the  center  bay  denoted  by  the  coordinate  Wj2  In  Figure  1  ->d  the 
supporting  structure  surrounding  ^fhe  center  bay.  The  dynamic  characteristics  of  such 
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(b)  SUBSTRUCTURE  GEOMETRY  AND  NOMENCLATURE 
FIGURE  1 .  NINE-BAY  FLAT  STIFFENED  PANEL 


structure  are  extremely  comolex.  Analytically,  several  methods  of  analysis  have  been 
applied  to  such  structure.  The  finite  element  approach  has  attracted  much  attention 
(References  3,  4,  and  5);  however,  this  approach  is  not  very  suitable  for  rapid  design 
application,  since  extensive  computer  facilities  are  required.  Clarkson  has  pointed  out 
that  an  elementary  approach  can  provide  estimates  at  least  within  the  scatter  of  fatigue 
data  (Reference  1) . 


1  .  Supporting  Structure 

The  common  stiffening  members  used  in  aircraft  construction  can  be  broadly  classified  as 
thin-walled  open-section  beams.  The  difference  between  this  type  of  structure  and  the 
elementary  beam  usually  encountered  in  practice  is  that  warping  of  the  cross-section 
normal  to  its  plane  must  be  considered  in  addition  to  the  axial  strains  introduced  in  the 
elementary  bending  theory.  The  discussion  here  shall  focus  on  the  torsional  restraint 
offered  by  such  stiffeners,  since  this  is  the  basic  mode  of  deformation.  Strictly  speaking, 
the  bending  and  torsion  motion  of  such  structure  are  coupled,  but  for  the  stiffeners  con¬ 
sidered  here,  the  bending  stiffness  is  several  orders  of  magnitude  greater  than  the  torsional 
rigidity.  Vlasov,  Reference  6,  and  Oden,  Reference  7,  give  very  complete  developments 
for  the  analysis  of  thin-walled  open-section  beams. 

Figure  2  illustrates  the  stiffener  geometry  under  consideration .  The  attach  point  is  taken 
as  the  point  at  which  the  stiffener  is  connected  to  the  panel.  It  is  assumed  that  the  stif¬ 
fener  is  free  only  to  rotate  about  the  attach  point  so  thar  it  is  necessary  to  present 
expressions  for  the  warping  constant  and  polar  moment  of  inertia  referenced  to  the  axis 
system  at  the  attach  point.  It  can  be  shown  that  (References  6  and  7)  the  effective 
warping  constantsT*  and  T*  are  given  by  the  expressions  (see  Figure  2) 

x  y 

r*  =r  +  s2i  -  2S  s  i  +  s2i 

x  e  z  zz  y  z  yz  x  yy 

for  stiffeners  parallel  to  the  x-axis  and 

r*=r  +  s2i  -  2S  s  i  +s2i 

y  e  z  zz  x  z  xz  x  xx 

for  stiffeners  parallel  to  the  y-axis. 

The  polar  moment  of  inertia,  referenced  to  rotations  about  the  attach  point,  are 


=  1  +1  +<c2+c2)a 

px  yy  zz  y  z 

(8a) 

=  1  +1  +(C2  +  C2)A 

py  xx  zz  x  z 

(8b) 

These  expressions  will  be  used  subsequently  for  developing  expressions  for  the  potential 
and  the  kinetic  energy  for  the  frequency  analysis  of  the  stiffened  panel. 


(7a) 


(7b) 
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2.  Frequency  Analysis 

The  development  here  assumes  that  the  cover  sheet  is  supported  at  the  stiffener  attach  line 
so  that  the  only  motion  of  the  structure  consists  of  transverse  displacement  of  each  panel 
bay  and  rotation  about  the  stiffener  attach  line.  Adjacent  panel  bays  are  constrained  so 
that  slope  continuity  across  the  stiffener  attach  line  is  preserved.  Energy  methods  are  used 
to  obtain  expressions  for  the  modal  stiffness,  modal  mass,  and  the  generalized  force. 

2.1  Cover  Sheet  -  For  the  rectangular  panel  illustrated  in  Figure  1,  the  transverse  dis¬ 
placement  function  for  a  bay  with  dimensions  a;  and  bj  is  denoted  by  w;j  (x,y)  where 
(x,y)  denotes  a  local  axis  system  for  each  panel  bay.  The  cover  sheet  is  assumed  to  be 
of  constant  thickness  and  uniform  material  properties. 


Assuming  a  mode  for  each  panel  bay  of  the  form 


w  (x,y)=C  sin(=^)Sm(^) 


(9) 


it  is  seen  that  there  are  nine  amplitude  parameters,  C..,  to  be  determined .  For  bay  motion 
described  by  Equation  9  there  are  eight  independent 'slope  continuity  conditions  across 
the  stiffener  attach  line  between  adjacent  bays.  Selecting  the  center  bay  amplitude  as 
the  reference  amplitude,  the  parameters,  C.. ,  given  in  Equation  9  are  expressed  in  terms 
of  the  center  bay  amplitude  as  indicated  in''  Table  I .  The  center  bay  motion  is  given  as 


(10) 


TABLE  I 

VALUES  OF  C../-W  FOR  A  NINE  BAY  STIFFENED  PANEL 


i  =  l 

2 

3 

i=  i 

(-l^Hfojbj/a^ 

(-i)n(b1A2) 

(-Om^ajbj/ajbj) 

i  =  2 

(-1  )m(a,/a2) 

1 

(-l)m(a1/a2) 

i  =  3 

(-l)m+n(a1b1/a2b2) 

(-l)"(b|A2) 

(-l)m+n(a  jbj/ajbj) 

For  a  thin  rectangular  plate  of  dimensions  a  x  b  undergoing  bending  deformations,  the 
expression  for  the  kinetic  energy  is 

a  b 

T  /  J  w2(x,y)dydx  (11) 

o  o 
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on  ’  he  potentiol  energy  is 

a  b 

U  =  -Id  f  f  [w2  +  w2  +  2vw  w  +  2(1  -  v)w2  1  dydx  (12) 

2  J  J  L  ,xx  / yy  ,xx  ,yy  ,xyj 

0  0 

For  a  simply  supported  plate  with  dimensions  a  x  b  vibrating  in  the  (m,n)^  mode,  the 
expressions  for  the  kinetic  and  potential  energy  become 


T  =  lYtabW2 
o  mn 

(13) 

, .  _  1  n4D  2  2.  ,,  .  w2 

U-87Tm"  fmn(b'°,Wmn 
where  2 

'«.«>••> -[(SXs)4©©] 

(14) 

and  the  mode  shape  is  given  by 

w(x,y)»sin(2H)sin(^)Wmn 

(15) 

Substituting  the  mode  functions  given  by  Equation  9  with  amplitudes  given  in  Table  1  into 
Equations  13  and  14  and  summing  the  terms  from  each  bay,  one  obtains  the  expression 
for  the  kinetic  energy  of  the  cover  sheet  as 

*1)2*1 

(16) 

and  the  potential  energy  as 


Up  *55^m2n2  tFmn(b2'a2l  +  2(blA2)  +  2(ol/o2)  F1Jb2'al) 

+  4(°)/o2>(b|'/b,2)  ^1'°,))  2WL 
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Equations  16 and  17  will  be  used  to  formulate  the  modal  mass  and  stiffness  of  the  structure. 

2.2  Supporting  Structure  -  The  substructure  geometry  and  nomenclature  are  illustrated  in 
Figure  1 .  Figure  2  illustrates  the  cross-section  geometry  for  the  stiffeners  indicating  the 
relative  location  of  the  shear  center,  attach  point,  and  centroid  for  the  stringers  parallel 
to  both  the  x-axis  and  the  y-axis.  The  stiffener  rotation  is  defineJ  as 
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,  .  .  ( m  TTX\  _ 

9..(x)  =  sin  ( - )  ©.. 

il  \  a.  /  i| 

1  I 

9n(y)  =  sin  (t*)  @ii 


where  i  is  the  stiffener  number  and  j  is  the  bay  index  number. 

Assuming  that  each  stiffener  is  uniform  along  its  length,  the  kinetic  energy  for  stiffeners 
parallel  to  the  x-axis  is 

3  a. 

T  .  =  i  (v  I  ) -  Y  I  1  e2.(x)dx  i  =  3,4  (19 

si  2  px  i  L  J  i  j 

i=i  o 

and  for  stiffeners  parallel  to  the  y-axis 


3  b. 


T.=i(yl  )j7  I  e..(y)  dy 
si  2  py  i  4_i  J  1 1 


=  1/2 


where  I  and  I  are  given  by  Equation  8  . 
px  py 

The  potential  eneigy  of  stiffeners  parallel  to  the  x-axis  is  given  by 

3  a.  a. 

Vll  <“?.  /V.W( GJX>.  /  ejfwdx 

i=i  o  11  o 

and  for  the  stiffeners  parallel  to  the  y-axis 

3  b.  b. 

usi4l  <E7i  / '  + (GVi  /  9l^(y)dy 


where  T*  and  V*  are  given  by  Equation  7  . 
x  y 


The  stringer  twist  amplitudes,  0..,  normalized  to  the  center  bay  plate  amplitude,  ^ 
are  given  in  Table  II. 
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TABLE  II 

VALUES  OF  0..AW 

i  j  2  mn 


i*  ^ 

2 

3 

i  =  1 

(-l)nmTTb1/a2b2 

nrm/a2 

(-l)nmtTb1/a2b2 

i  =  2 

(-1)m+nmTTb^/a2b2 

(-l)mmn/a2 

(-l)m+nrr,nb^/a2b2 

i  =  3 

(-1)mnna1/a2b2 

CM 

< 

1= 

C 

(-l)mnna1/a2b2 

i  =4 

(-l)m+nn  TTOj/c^bj 

(-1)nnn/b2 

(-l)m+nn  TTd  ^0^2 

Substituting  Equation  18a  into  Equations  19a  and  20a  and  Equation  18b  into  Equa¬ 
tions  19b  and  20b  ,  performing  the  integrations,  and  substituting  the  amplitude 
coefficients  given  in  Table  II,  the  expressions  for  the  kinetic  and  potential  energy  become 


i  =  3,4  (21a) 


U  .  = 

si 


6  (Er*). 

TT  X  I 

4  "  375" 

°2b2 


4  2 
m  n 


for  stringers  parallel  to  the  x-axis,  and 


T  . 

si 


i  =  3,4  (21b) 


i  =  1,2  (22a) 


u.=x 

si  4 


<EI7i  2 

IT" 


°2b2 


!  4  r 

n  L 


1  +  k  .9  +  2 

y  * 2 


for  stringers  parallel  to  the  y-axis 

where  k2..  =  (GJ  /Er*).  (b./nn)2;  k2..  =  (GJ  /Er*).  (a./mn)2. 
yi|  /  y  i  |  XI I  x'  x  |  f 


1  =  1,2  (22b) 


For  simplicity,  it  will  be  assumed  that  the  stiffeners  parallel  to  the  x-axis  are  identical 
and  that  the  stiffeners  pcrallel  to  the  y-axis  are  identical.  Then  Equations  21a  and  21b 
become 
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4 


2  2  a, 

(y  , 

px  \j 2  '”2' 

6  (Er*) 


ii  tt'  x'  4  2  r  ,2  0 

u»  =  T1SJm"  Ll  +  kx2+2 

°2b2 
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and  Equations  22a  and  22b  become 

2  2  b 


bl  3 


Vt*»V^[,  +  2(i5)] 

6  (Er*) 


"  '  tm2n4ri+k2,+  2 


U  =  ^ - *-4r  m  n  £ 


sy 


4  “O' 
°2b2 


y2 


fe( 


2w2 

(23a) 

2  mn 

!('  +  k«l)]  2Wmn 

(23b) 

,w2 

(24a) 

2  mn 

l('  +  kyl)]  2Wi 

(24b) 

where  k2.  =  (GJ/Er*)x  (a./mn)2;  k2.  =  (GJ/Er*)y  (b./ntr)2. 

2.3  Modal  Stiffness  and  Modal  Mass  -  With  reference  to  Equations  4  and  5,  the  modal 
stiffness  and  modal  mass  for  the  (mfn)^  mode  of  the  structure  are  obtained  by  summing 
the  contributions  of  each  component.  From  Equations  16  and  17  and  Equations  21  and 
22  the  expression  for  the  modal  stiffness  is 


The  expression  for  the  modal  mass  is 
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2  o  3 

Mmn  =  ?{y,(°2  +  2’l)(b2  +  Sk))  +  'px*  (bf>  ['  +  2  (^)  ] 


V' 


+  -2<y'p^(^['  +  2^)3]} 


(26) 


The  frequency  for  the  (m,n)^  mode  is  given  by 


f  =  J- V(K  /M  ) 
mn  2rr  mn  mn 

2.4  Generalized  Force  -  For  a  time-dependent  pressure  distributed  over  the  surface  of  the 
panel,  the  expression  for  the  generalized  force  is 


Pmn(f)  =  /  P(*,y,t)  w(x,y)  dA 
A 


(27) 


Assuming  that  the  pressure  is  uniformly  distributed  over  the  surface  of  the  structure,  sub¬ 
stituting  the  mode  functions  for  each  bay  and  integrating  over  the  surface  of  the  structure, 
the  expression  for  the  generalized  force  is 


P  (l)  =^[(-!)m  -  1][(-Dn  -  OR,  P« 

mn  L  mn  o 

rr  mn 


(28) 


where  R„„  -  1  *  2(-)>"  ♦  2(-.f  ♦  4(-l)m+n 


3.  Stiffener  Response  to  Uniform  Static  Pressure 

For  a  simply  supported  panel  loaded  by  a  uniform  pressure,  q  ,  the  pressure  can  be 
resolved  into  a  Fourier  series  as 


P  = 


16q, 


tt  mn 

For  a  displacement  function  of  the  form 


rv  o  .  /mrrx\  .  (n tt y * 

LL~t~  \rri 


(29) 


(30) 


the  static  response  of  the  panel  to  the  load  qQ  in  the  (m,n)^  mode  is 
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(31) 


W 

mn 


2  2 

1 6q  a  b 
o 


6  3  3  _  _  (u  \ 
tt  m  n  D  F  (b,a) 
mn 


The  stress  response  to  the  static  loading  must  consider  stresses  in  the  cover  sheet  and 
stresses  in  the  supporting  ribs.  Ballentine's  design  method  for  establishing  skin  thicknesses 
of  acoustically  loaded  panels  (Figure  73,  Reference  8)  is  sufficiently  accurate  so  that 
additional  effort  on  this  aspect  of  the  problem  is  not  required  for  stiffened  panels. 

The  stress  response  of  the  substructure  is  of  primary  interest.  The  static  loading  of  the 
panel  is  assumed  to  be  transferred  to  the  stiffener  by  a  shear  force.  The  intensity  of  the 
load  transmitted  can  be  estimated  from  thin  plate  theory  by  computing  the  shear  loading 
along  the  edge  of  the  panel.  Fora  panel  deforming  as  given  in  Equation  30,  the  distri¬ 
buted  shear  loading  along  the  edge  x  =  0  is  given  by  the  expression 


16qob[Fln2(b'a)  +  0  ‘v)GXe)]  .  /nrry\ 
v(°,y)  = - 3  5 - - - - - sm  V  b  j 


tt  mn  F  (b,a) 
mn 


(32) 


for  a  uniform  loading  q  .  The  resultant  of  this  shear  loading  is  obtained  by  integrating 
Equation  (32)  over  the  °length  of  the  edge  and  subtracting  the  two  corner  reactions  o(  the 
plate  (Reference  9,  pp.  105,107).  The  resultant  shear  load  is  given  as 


16qob2[(-l)n-l] 
4  171727-— 

tt  mn  F  '  (b,a) 


(33) 


Considering  the  stringer  to  be  clamped  on  both  ends  (for  bending)  and  loaded  by  a  uniform 
load  of  intensity  V/b,  the  bending  moment  at  the  support  has  magnitude  (setting  qQ  =  1) 


M  4b3  [(-1)"  -  I] 

o  7~4  3.1/2,.  . 

3tt  mn  F  (b,a) 


(34) 


The  bending  stress  in  the  flange  of  the  stiffener  is  given  in  Reference  7,  p.  193,  as 


a 

o 


I  M  h 
zz  o 


I 

XX  zz 


M  h 
0 

21* 


(35) 


where  h  =  height  of  the  stiffener  cross-section. 

The  area  moments  of  inertia  given  in  Equation  35  are  referenced  to  the  centroidal  axes. 
Expressions  for  the  area  moments,  St.  Venant's  torsion  constant,  and  the  warping  constant 
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are  given  in  Appendix  I  for  constant  thickness  zee,  channel,  and  hat  cross-section 
shapes. 


D.  Box  Structure 

The  broad  definition  of  a  box  type  structure  is  taken  to  be  two  parallel  plates  stiffened  by 
deep  shear  ribs.  Such  structure  can  be  taken  as  a  model  for  flap,  fin,  or  pylon  type  struc¬ 
ture  encountered  on  modem  aircraft.  The  box  structure  shall  be  classified  by  the  number 
>  of  cells  formed  by  the  intersecting  cover  sheets  and  ribs.  Two  configurations  shall  be 

considered  for  the  frequency  analysis:  three  cell  and  nine  cell  box  structures.  These 
structures  are  illustrated  in  Figure  3.  The  three  cell  configuration  has  been  investigated 
by  Clarkson  and  Abrahamson,  Reference  10,  and  by  Sen  Gupta  and  Mead,  Reference  11  . 

The  presentation  here  follows  that  of  Reference  10  although  the  development  was  independent. 


1  .  Derivation  of  Mode  Shapes 


The  detailed  derivation  of  the  mode  shapes  for  the  Rayleigh  frequency  analysis  of  the 
three  cell  box  structure  is  presented  in  this  section.  Consider  the  structure  illustrated  in 
Figure  3a  with  the  coordinate  nomenclature  for  each  component.  The  mode  functions  for 
each  component  are  taken  as  (local  x-y  coordinates  are  implied  for  each  component): 


jWj(x,y)  = 

=  sin 

(^) 

sin 

(cSfi 

W 
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jw2(x,y)  = 

=  sin 

(¥) 

sin 
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Vj(x,z)  = 
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(Vs) 
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(v1) 

V. 

v2(x,z)  = 

=  sin 

(v=) 

sin 

V2 

2w1(*/y) s 

=  sin  1 

(Vs) 

sin 

(?Sf) 

2Wt 

2w2(x,y) s 

=  sin  1 

(Vs) 

sin 

(t?) 

w 

2  2 

2w3(*,y) 8 

:  sin  | 

(Vs) 

sin 

W 

2  3 

(36) 


where  i,  k,  K>,  m,  n,  p,  q,  r,  and  s  are  integer  mode  numbers. 


(a)  THREE  CELL  CONFIGURATION 


(b)  NINE  CELL  CONFIGURATION  (EXPLODED) 
FIGURE  3.  BOX  STRUCTURE  CONFIGURATIONS 


/ 


/ 


Assuming  that  each  joint  is  allowed  only  to  rotate  and  that  all  components  intersecting  at  a 
joint  hove  the  same  slope  (zero  shear  condition),  one  obtains  the  constraint  conditions: 


,w0 
1  2/y 


,U>l)  +  Vj  2<0)  =  0 

lw2,y(bl)  +  v2,z(0)=0 

(0)  +  V.  (0)  =  0 
*  /  * 

+  “2,z<°>=0 

cr 

+ 

< 

N 

_3~_ 

II 

o 

2w2,y(bl)  +  v2,z<h)=° 

(0)  +  v.  (h)  =  0 

1  1  ^ 

2W3,y(0>  +  V2,z:(h)  =0 

(37) 


The  eight  homogeneous  equations  given  above  have  a  non-zero  solution  (i.e.,  the  ampli¬ 
tudes  jWj,  etc  . ,  are  non-zero)  only  if  the  determinant  of  the  coefficient  matrix 

vanishes.  Substituting  from  Equation  36  into  Equations  37  and  evaluating  the  determinant 
of  the  coefficient  matrix,  one  obtains  the  condition  on  the  mode  numbers  that 
(-l)n+r  =  (-iyt+p  for  the  determinant  to  vanish.  That  is,  the  mode  numbers,  n,  r,  l,  and 
p  must  be  selected  such  that  (n+r)  and  (6+p)  are  either  both  odd  or  both  even  (n  =  1 ,  r  =  1 , 
1=2,  p=1  would  not  be  a  valid  mode  number  set,  for  example). 


The  above  derivation  is  too  general  to  be  of  design  value  and  has  been  presented  only  to 
illustrate  the  general  technique  for  establishing  valid  modes  for  the  Rayleigh  energy 
analysis. 

2.  Three  Cell  Box  Structure 


The  derivation  of  expressions  for  the  modal  stiffness  and  modal  mass  is  rather  straight¬ 
forward  once  the  mode  shapes  for  the  structure  are  established .  For  design  purposes,  it 
will  be  assumed  that  the  rib  spacing  is  uniform  (b^  =  =  b^  =  b  in  Figure  3a),  that  the 

material  properties  of  the  cover  sheets  and  ribs  are  identical  and  that  the  only  modes 
considered  are  described  by  the  set  (p,m,r)  for  the  (x,y,z)  coordinate  directions. 

For  the  above  assumptions,  one  makes  the  following  substitution  for  mode  number  s  as 
given  in  Equation  36 


i  =  p;  k=6  =  q  =  r  =  s  =  m;  n  =  p  =  r  . 

With  this  assignment  the  determinant  of  the  coefficient  matrix  vanishes  identically. 
Expressing  the  eigenvector  of  component  amplitudes  in  terms  of  the  amplitude  ^  (center 
bayof  the  cover  sheet  in  the  plane  (x,y,0)),  one  obtains: 

1W1  =  H)f"  1W2  2W1  =  H)m+r  1W2 

1W3  =  Hr  ,W2  2W2  =  (-l)r  ,W2  (38) 

Vl  =  '  (b)  (?)  1W2  2W3  =  (-l)m+r  ,W2 
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v  <-'r+1  (£)(?)  ,w2 


(38  cont'd) 


Figure  4  illustrates  six  possible  mode  shapes  for  the  structure  for  the  indicated  modes  for 
a  section  through  the  structure  in  the  plane  (a/2,y,z).  Figures  4a,  4b,  4c  illustrate  the 
fact  that  for  fundamental  modes  in  the  cover  sheet  (p  =  m  =  1)  all  rib  modes  (r  =  1,2,3. . .) 
can  be  excited.  Also,  if  r  is  odd,  the  two  cover  sheets  are  out-of-phase  and, if  r  is  even, 
the  cover  sheets  are  in-phase.  Since  the  excitation  pressures  are  on  the  cover  sheets  and 
the  generalized  force  will  be  maximum  for  the  fundamental  mode  in  the  cover  sheet,  it 
should  be  possible  to  excite  any  rib  mode.  This  effect  has  been  observed  experimentally. 

2 . 1  Frequency  Ana  lysis  -  The  kinetic  energy  and  potential  energy  expressions  are  obtained 
by  substituting  the  mode  shapes  for  each  component  (as  defined  by  Equations  36  and  38) 
into  Equations  11  and  12.  The  expression  for  the  modal  stiffness  is 


Kpm,=  £  VV{Va>)  +  (cXr/  Frp<0'h)} 


and  for  the  modal  mass 


Mpn,r'ps°b{3  +  (&>  (f)(7)  } 


where  it  has  been  assumed  that  bj  =  bj  -  b^  =  b. 

The  frequency  of  the  (p,m,r)^  mode  is 

f  K  /M  ),  Hz 

pmr  2tt  pmr  pmr 

2.2  Generalized  Force  -  The  expression  for  the  generalized  force  for  the  three-cell 
box  structure  will  be  derived  assuming  uniform  pressures  of  different  intensities  over  the 
two  cover  sheets.  Denoting  the  pressure  on  the  upper  cover  sheet,  (x,y,0)  plane,  as 
Pj(t)  and  the  pressures  on  the  lower  cover  sheet,  (x,y,h)  plane,  as  P2(t)  and  assuming 
that  the  pressure  is  positive  if  directed  toward  the  interior  of  the  structure,  the  expres¬ 
sion  for  the  generalized  force  is 

3  a  b  3  a  b 

Ppmr^  =  Pl^Z  /  /  iw|(x^)dydx  +  P2W  Ya  /  /  2wi^x,y^dyclx  (4°) 


i=1  0  0 


i=1  0  0 


Substituting  for  the  modes  shapes  jW.(x,y)  and  2w-(x#y)  and  noting  the  relations  given  by 
Equations  36  and  38,  one  obtains 

p  (t)  =-v-[p,(t)  -(-i)rp,(t)Tl  -  2H)mlH)m  -  |T(-1)P  - 1]  <«1> 

h  TT  run  J 
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FIGURE  4 


MODE  (1,1,1) 


MODE  (1,3,2) 


MODE  SHAPES  FOR  A  THREE  CELL  BOX  STRUCTURE 
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Equation  41  clearly  illustrates  that  only  modes  for  p  and  m  odd  will  be  excited,  as  would 
be  expected . 

3.  Nine  Cell  Box  Structure 

The  configuration  of  the  nine  cell  box  structure  is  illustrated  in  Figure  3b.  This  configu¬ 
ration  was  adopted  for  the  experimental  program  described  subsequently.  The  displacement 
functions  for  each  component  are  denoted  in  Figure  3b  where  cover  sheet  displacement  is 
given  by  the  functions  ^w..  and  2wj.  (positive  in  the  positive  z  direction)  and  rib  displace¬ 
ment  functions  jU.,  and  jvj/  2vj  (positive  in  the  positive  x  and  y  directions,  respec¬ 
tively)  .  As  described  for  the  three  cell  box  structure,  the  application  of  a  zero  shear 
condition  at  each  joint  of  a  rib  and  cover  sheet  allows  one  to  describe  a  mode  in  terms  of 
a  mode  number  set  and  a  single  amplitude.  Also,  the  analysis  presented  here  assumes  that 
a  mode  is  described  by  the  mode  numbers  (m,n,q)  so  that  the  assumed  modes  for  each 
component  have  the  form 

kw  <x,y)  =  sin  (=£*)  sin  (=£*)  fcW 
'  I* 

kU;(y,z)  =  sin  (^)  sin  kU.  (42) 

i 

kvi(x,2)  =  sin  (=£2)  sin  (Sjji)  kV. 

I 

for  the  component  of  the  cover  sheet  or  rib.  Selecting  the  center  bay  of  the  upper 
cover  sheet  as  the  arbitrary  amplitude,  -^22'  amPl'fudes  °f  the  components  in  terms 
of  the  mode  number  and  geometry  are  given  as: 

lW12  ‘  1 W31  '  lWi3  '  iW33  ‘  <V°2><b,A2>  l\,n, 

<«■> 

,W,2=lW32  =  <-'>mV°2>  1W™, 

for  the  upper  cover  sheet  amplitudes; 

,V,u3  =  (h/ojKyv  ^ 

,U2=  -  (h/ojXm/q)  lWmnq  (43b) 

2Ui  Ui 


i 


(43c) 


lV,-lV3-(“,)m^1^  1Wmnq 
,V2  =  -  (h/b^n/q)  jWmnq 


2Vi  =  <-'>  ,Vi 
for  the  rib  amplitudes;  and 

w  =  W 
2  ij  1  ;  1  Wij 

for  the  lower  cover  sheet  amplitudes  and  =  jW22 


(43d) 


3 . 1  Frequency  Ana  lysis  -  For  the  frequency  analysis,  it  will  be  assumed  that  the  material 
properties  of  all  components  are  identical  and  that  the  thickness  of  the  components  are  as 
illustrated  in  Figure  3b.  The  modal  stiffness  and  modal  mass  for  the  (m,n,q)^  mode  are 
obtained  by  substituting  the  component  mode  shape  as  given  in  Equation  42  with  amplitudes 
given  by  Equations  43  into  Equations  11  and  12  and  summing  for  all  components.  Neglecting 
the  detailed  algebra,  the  expression  for  the  modal  mass  is  given  by 


M 

mnq 


(44) 


where  .  p  =  y ,  t 
Is  Is 


=  1  +  2(a1/a2)3  +  2(b1/b2)3  +  (bj/b^' 

V,  +  aV“:H7)(fe2[i  +  *VVl 

1  r  2 


The  expression  for  the  modal  stiffness  is  given  as 

4 

-  ,D  0  0  ,t  3  ,t  3 

2xf m  n  ([' + (7;  ] K, +  2  (7) 


mnq  a 


T 


where  =  E  |t3/12(l  -v2) 


Ks  =  Fmn(b2'°2)  +  2[(al/°2)Fmn(b2'0l)  +  (bA)Fmn(bl 
+  2(a1/a2)(b1/b2)Fmn(b1,a1)] 


(45) 
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Kf  =  (h/bj)  [Fmq(h,02)  +  2(a1/a2Fmq(h,a1)] 

+  »’/«2><2,/lV3  CF„q(h-b2)  +  2(bl/b2)Fnq<h'bl):1 


The  frequency  of  the  (m,n,q) 


th 


mode  is  given  by  the  expression 


f 

mnq 


J-V(K  /M  ),  Hz 
2tt  mnq  mnq 


3.2  Generalized  Force  -  The  relationship  for  the  generalized  force  for  the  nine-cell  box 
structure  is  obtained  in  a  manner  identical  to  that  described  for  the  three  cell  box  struc¬ 
ture.  Assuming  that  a  pressure  is  positive  if  it  is  diiected  towards  the  interior  of  the 
structure  and  denoting  the  pressure  on  the  upper  surface  as  p^(t)  ((x,y,C)- plane)  and  on  the 
lower  surface  as  p2(t)  ((x,y,h)- plane),  the  generalized  force  for  the  (m,n,q)^  mode  is 

a2b2  a 

P(t)  =  -r1  <*<*> "  H>q  <46> 

mnq  z  I  z  mn 

TT  mn 

where  T  =  [(-i)™  -  1][(-l)n  -  1]  R 

mn  mn 

R  is  defined  in  Equation  28. 
mn  n 

Equation  46  clearly  illustrates  that  only  modes  for  m  and  n  odd  will  be  excited,  as  would 
be  expected . 


4.  Stress  Response 


As  indicated  by  Equation  3,  the  stress  response  of  the  skin  and  ribs  for  box  structure  due  to 
a  uniform  static  pressure  of  unit  magnitude  is  required  to  estimate  the  mean  square  or 
R.M.S.  stress  response  of  the  structure.  The  approach  taken  here  is  dependent  upon  the 
flexibility  of  box-type  structure  as  described  in  Reference  1  and  is  essentially  different 
from  the  approach  taken  by  Ballentine,  Reference  8.  The  approximate  method  for  com¬ 
puting  stresses  in  continuous  slabs  as  described  by  Timoshenko  (Reference  9,  pp.  236-245) 
is  presented  here  —  modified  to  account  for  the  rib  supporting  the  cover  sheet.  This 
modification  was  first  described  by  Clorkson,  Reference  1 . 

Timonshenko's  approximate  method  relies  upon  the  synthesis  of  elementary  solutions  for 
uniformly  loaded  thin  rectangular  plates  to  obtain  estimates  of  stresses  in  continuous  slabs. 
The  basic  procedure  is  to  divide  the  bays  of  the  structure  into  components  such  that 
boundaries  between  adjacent  bays  are  considered  to  be  clamped.  Considering  the  plan- 
form  of  the  cover  sheet  for  three  cell  and  nine>cell  box  structure  as  illustrated  in  Figure  3, 
the  bays  of  structure  are  taken  with  the  boundary  conditions  as  illustrated  in  Figure  5. 
Timoshenko  presents  tabulated  data  for  estimating  the  bending  moments  at  the  panel  center 
and  the  center  of  a  clamped  edge  for  uniformly  loaded  plates  with  various  boundary  con¬ 
ditions.  At  a  point  on  the  boundary  between  two  plates,  the  resultant  bending  moment  is 
approximated  by  the  average  value  of  the  component  bending  moments.  As  illustrated  in 
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a)  BOUNDARY  CONDITIONS  FOR  THREE  CELL  BOX  STRUCTURE 


b)  BOUNDARY  CONDITIONS  FOR  NINE  CELL  BOX  STRUCTURE 


-  SUPPORTED  EDGE 

CLAMPED  EDGE 


FIGURE  5  COMPONENT  BOUNDARY  CONDITIONS:  BOX  STRUCTURE 
STRESS  CALCULATION 


22 


Figure  5,  four  different  sets  of  boundary  conditions  are  to  be  considered.  Case  1  repre¬ 
sents  a  rectangular  plate  with  three  edges  supported  and  one  edge  clamped .  Case  2 
represents  a  rectangular  plate  with  two  opposite  edges  supported  and  two  opposite  edges 
clamped .  Case  3  represents  a  rectangular  plate  with  one  edge  supported  and  three  edges 
clamped.  Finally,  Case  4  represents  a  rectangular  plate  with  all  edges  clamped.  The 
bending  moments  at  the  center  of  a  clamped  edge  for  a  rectangular  plate  with  dimensions 
a  x  b  uniformly  loaded  by  a  pressure  of  magnitude  qQ  is  given  by  Timoshenko  as: 


Case  1 

M  =  60  q 
y  2  no 

at  (x,y)  =  (a/2,0) 

Case  2 

M  -  60  q  -t2 
y  3  o 

at  (x,y)  =  (a/2,0)  &  (a/2,b) 

Case  3 

M  =  y  q 
x  5  o 

at  (x,y)  =  (0,b/2)  &  (a , b/2) 

M  =  6C  q  l? 
y  5  a 

at  (x,y)  =  (o/2,0) 

Case  4 

M  =  y.  q  -L2 
x  6  no 

at  (x,y)  =  (0,b/2)  &  (a, b/2) 

M  =  6 .  q  -L2 
y  6  no 

at  (x, y )  =  (a/2,0)  &  (a/2,b) 

where  L  is  the  smaller  of  the  dimensions  a  and  b  and  the  subscripts  on  the  y's  and  6's  (the 
bending  moment  factors)  follow  Timoshenko's  notation.  Plotted  values  of  the  bending 
moment  factors  versus  panel  aspect  ratio  b/a  are  present  in  Figures  6,  7,  and  8.  For  more 
accurate  values  of  the  bending  moment  factor,  the  referenced  tabulated  values  can  be 
used.  The  asymptotic  values  for  the  bending  moment  factors  are  listed  in  Table  III . 


TABLE  III 


ASYMPTOTIC  VALUES  OF  THE  BENDING  MOMENT  FACTORS 


b/a 

62 

63 

y5 

65 

y6 

66 

0.0 

.1250 

.0833 

- 

- 

.0571 

.0833 

00 

.1250 

.1250 

.0833 

.0566 

.0833 

.0571 

NOTE:  All  values  for  the  bending  moment  factors  as  presented  here  are 
positive  instead  of  negative  as  presented  by  Timoshenko. 


The  use  of  the  plotted  or  referenced  tabulated  data  for  computing  bending  moments  in  the 
cover  sheet  or  rib  of  the  panel  arrays  shown  in  Figure  5  is  best  illustrated  by  an  example. 


Example:  (a)  Three  Cell  Box  Structure  -  Considering  the  panel  array  illustrated  in  Figure 
5a,  the  values  for  and  6^  are  obtained  for  the  particular  value  of  panel 
aspect  ratio,  b/a  (Figure  6)  and  the  expression  for  the  bending  moment  at 
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l  = 

SMALLER  OF 

SPANS  o  AND  b 

b9 


point  A  (Figure  5o)  is  simply 


M  =(60  +  6Jq  (47) 

y  c  o  o 

where  -t,  is  the  smaller  of  the  two  dimensions  a  and  b 
N  =  2  for  a  stiffened  panel  array 
N  =  3  for  a  box  type  structure.  Ref.  1  . 


(b)  Nine  Cell  Box  Structure  -  Considering  the  panel  array  illustrated  in  Figure 
5b,  the  value  of  6^  is  obtained  for  the  two  panels  with  dimensions  02  x  b^ 
and  a,  x  b_  and  the  values  of  y  .and  6,  are  obtained  for  the  center  bay  so 
that  for  point  A 


M  =  M?+6X?)q  /N  (48a) 

and  for  point  B  7  51  6  7  0 

k-hV'&W*  (48b) 


where  is  the  smaller  of  the  two  dimensions  02  and  bj 
's  smaller  of  the  two  dimensions  02  and  b2 

is  the  smaller  of  the  two  dimensions  a^  and  b2 


Stress  is  obtained  from  Equations  47  and  48  by  use  of  the  simple  expression 

a  =  6M/t2  (49) 

where  t  is  the  thickness  of  the  cover  sheet  or  rib. 


E,  Frequency  Analysis  of  a  Three  Cell  Wedge  Structure 
1 .  Basic  Considerations 


Wedge  structure  is  typical  of  aircraft  control  surfaces  such  as  ailerons,  trailing  edge  of 
flaps,  elevators  and  rudders.  The  presentation  here  is  a  very  simplified  approach  to  a 
complicated  structural  dynamics  problem.  The  objective  will  be  to  obtain  an  estimate  of 
the  natural  frequencies  of  such  structure.  Due  to  the  very  simple  mode  shape  assumed  for 
the  rib  in  the  radial  direction,  any  attempt  to  obtain  stress  estimates  would  be  of  dubious 
value.  However,  the  frequency  estimates  do  seem  reasonable. 

2.  Frequency  Analysis 

Consider  the  three  cell  wedge  structure  illustrated  in  Figure  9.  For  structural  modes 
symmetric  about  the  x-y  plane,  mode  shapes  for  a  component  of  the  cover  sheet  are 
assumed  to  be  of  the  form 
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/ 


w  (Xi.yl-sin^XjJiin^jWi 


vk(r,8) 


(50a) 

(50b) 


where  i  =  1,2;  j  =  1,2,3;  k  =  1,2, 3, 4;  q  =  1,3,5, ...  The  mode  functions  w..  represent 
displacements  normal  to  the  x;-y  plane  and  v^  represents  the  rib  mode  functions.  For  a 
zero  shear  condition  at  a  joint  of  a  rib  and  a  cover  sheet,  one  obtains  the  relations 


l^k 

r  90 


+ 

8= a 


^lj 

9y 


=  0 

y=0,b. 


<3Vk  3w^ 

30  Q_  +  9y  _n  . 
0--a  '  y-0,b. 


(51) 


Substituting  Equations  50  into  51,  the  rib  amplitudes,  V^,  are  expressed  ir.  terms  of  the 
cover  sheet  amplitudes,  W..,  as 

i| 


V,=2aH)^2"-1>/2(^)w)2  =  V. 
V2  =  2aH>(^(l)(0W12  =  V4 


W- =  <-'»"£) 


11 


w_.  =  -  w,. 

21  1 1 


2 


w  =  w 

12  13 


where  the  amplitudes  have  been  normalized  to 

The  geometry  of  the  sector  plate  is  illustrated  in  Figure  10. 


(52) 


The  expression  for  the  potential  energy  of  a  sector  plate  (rib)  is 


+  2(l-v)^v0J  ]jrdrd0  (53) 

$  r 

Substituting  for  the  mode  shapes  given  in  Equation  50b,  the  expression  for  the  potential 
energy  becomes 
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TABLE  IV 

VALUES  OF  THE  INTEGRAL  I 


m 


The  kinetic  energy  of  a  sector  plate  is  given  by  the  expression 

a  a 

T  =  IYtr  J  J  v2(r'd)rdrdQ  (55) 

-a  0 

Substituting  for  the  mode  shapes  given  in  Equation  50b,  the  expression  for  the  kinetic 
energy  is 


The  expressions  for  the  kinetic  energy  and  the  potential  energy  of  the  flat  rectangular 
cover  sheets  are  obtained  by  substituting  Equation  50a  into  Equations  11  and  12.  For  the 
normalized  values  of  amplitude  given  in  Equations  52,  the  expressions  for  the  kinetic  and 
potential  energy  are  expressed  in  the  form  of  Equation  39. 

Assuming  that  the  cover  sheets  and  ribs  are  made  of  the  same  material  and  that  the  ribs  all 
have  the  same  thickness,  the  expressions  for  the  modal  stiffness  and  mowil  mass  are 
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Kmnq  =  Tsf"  1°S  {('  +  (fij)  )[FJb2'o)  +  (bf)  Fmn(bl  'o)  +  (bf^  Fm„(b3'o)] 


V  3  /2a\2  /  S 


2  2 


+ (7)  ©  + 15 + ^ + 4© +4[©  - ']  u  <57> 

1  s  n  2 


rj'v^iO + frl' +  (etJ + 0] + t  (i)  0)  L'  •  dhll 


2  '  2'J  I  s 


where  ,D  =  E  ,t3/12(l  -  v2) 
Is  Is 

S  =  2aa 


The  frequency  of  the  (m,n,q)  mode  is  given  as 


f  =  J-V(K  /M  )  ,  Hz. 
mnq  2n  mnq  mnq 


Typical  mode  shapes  for  the  rib  modes  (1,1,1),  (2,1, 1),  and  (3,1,1)  are  illustrated  in 
Figure  1 1 . 


F.  Approximate  Frequencies  of  Cylindrical  Panels 
1  .  Szechenyi's  Approximate  Formulae 

Szechenyi,  Reference  12,  presents  an  approximation  technique  for  estimating  the  natural 
frequencies  of  curved  panels.  Szechenyi's  result  has  been  checked  against  analytical 
and  experimental  data  presented  by  Petyt,  Reference  13,  with  good  agreement.  The 
following  formulae  are  taken  directly  from  Reference  12. 

For  either  simply  supported  or  clamped  boundaries  of  the  rectangular  cylindrical  shell  ^ 
illustrated  in  Figure  12,  Szechenyi  obtains  the  expression  for  the  frequency  of  the  (m,n) 
mode  as 


tk2  tk2  19r  1/2 

m  ,  n  ,  12Gl 


,  I  ID  f/i  *  .  1  4/  .  m  ,  n  ,  1^1 
fmn  2tt  V  M  L  ^m  kn)  ax  D  °y  D  f2R2j 

where  D  =  Et3/12(l  -v2) 

M  =  yt 

m  =  mode  number  in  the  longitudinal  direction 
n  =  mode  number  in  the  circumferential  direction 
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MODE  (3,1,1) 


a  =  longitudinal  panel  edge  stress 
a  -  'hoop'  stress 

y 

The  terms  k  ,k  ,  and  G  depend  upon  the  type  of  boundary  conditions  of  the  panel.  For 
simply  supported  edges: 


For  clamped  edges: 


k  =  (mn/a) 


k  =  (n  tt/1d) 

il 

G  =  (1  -vVV  +  k2)2 
m  m  n 


k 

m 


=  mrr/C  a 
m 


k  =nn/C  b 
n  n 


(60) 


(61a) 


values  of  G  are  given  in  Figure  13  and  values  of  £  and  £  as  a  function  of  edge  stress 
parameter  are  presented  in  Figure  14. 

For  zero  edge  stress. 


c  *  1/(1  +.  15/m) 
m 

C  w  1/(1  +  .15/n) 
n 


(61b) 


In  terms  of  an  internal  pressure, p,  the  hoop  stress  is  given  by 

=  pR/t  (62) 

2.  Plumblee's  Approximate  Formulae 

Plumblee,  Reference  8,  has  developed  expressions  for  the  frequency  ratio  and  bending 
stress  ratio  for  a  thin  cylindrical  panel  compared  to  a  flat  panel .  The  boundary  conditions 
are  considered  to  be  clamped  .  For  the  (m,n)^  panel  mode,  the  frequency  ratio  is  given 
as 


UJ 

c 

u> 


0.006  b3/tR 
A4  +  0.61A2  + 


(63) 


and  the  stress  ratio  for  a  point  on  the  middle  of  the  straight  edge  is  given  as 


'  SIMPLY  SUPPORTED  SHELL  AND 
CYLINDERS  WITH  SIMPLY 
•SUPPORTED  ENDS  FOR  ALL 
MODE  NUMBERS  m  AND  n 


FIGURE  13  G  FACTOR  FOR  m 


r  4/2*2  1-3/4  [  /  A  +  0-034_\S-  (64) 

o  f  0.006b  /t  R  1  +  0.453(  — - T'7/* 

1  Va4+9.62A+1/  1 

,  .  .  ospec,  ralio .  The  constants  in  the  above  equations  hove  been 

,here  A  =  b/a,  the  panel  aspect  rat.o 

letermined  empirically.  restrictions  must  be  placed  upon 

rhe  Rayleigh-Ritz  -.rictions  ore  cp  <  0.35, 


Ill  -  EXPERIMENTAL 


A.  Introduction 


The  objective  of  the  experimental  program  was  to  establish  more  accurate  methods  for  the 
design  of  acoustic  fatigue  resistant  internal  airframe  structure.  The  analytical  effort 
described  in  Section  II  was  used  as  a  guide  for  the  specimen  design  and  interpretation  of 
the  observed  experimental  results. 

To  accomplish  the  intent  of  the  investigation,  high-intensity,  random-noise  acoustic 
fatigue  tests  of  36  specimens  were  conducted.  The  specimen  configurations  consisted  of  12 
stiffened  panel  designs  and  6  box  specimen  designs  with  2  specimens  of  each  design  to 
check  repeatability.  The  parameter  range  covered  by  the  specimen  designs  was  as 
follows: 


(a)  Stiffened  panel  specimens: 

skin  thickness 
rib  thickness 
rib  spacing 
aspect  ratio 
rib  shape 

(b)  Box  specimens: 

skin  thickness 
rib  thickness 
rib  spacing 
height 

panel  aspect  ratio 

Due  to  the  limited  number  of  test  speci 
pletely  covered .  A  description  of  the 
given  in  the  following  sections. 


0.032  to  0.050  inch 
0.025  to  0.040  inch 
6.0  to  9.0  inches 
2.0  to  3.0 

zee,  channel,  and  hat  cross-section 

0.040  to  0.063  inch 
0.020  to  0.063  inch 
7.0  to  10.0  inches 
7.0  and  10.0  inches 
2.0  to  4.0 

is,  the  above  parameter  range  was  not  com- 
t  specimens,  test  procedure,  and  test  results  is 


B.  Test  Specimen  Design 

The  test  specimen  design  wasaccomplished  in  two  phases.  For  each  phase,  6  stiffened 
panel  and  3  box  specimen  designs  were  considered.  The  first  group  of  specimens  was 
designed  using  the  applicable  design  nomographs  of  Reference  8.  The  second  group  was 
designed  using  the  same  nomograms  and  data  from  the  first  phase  experimental  program. 
Fastener  edge  distances  and  other  design  details  were  selected  in  accordance  with  accepted 
aircraft  standards.  Standard  manufacturing  techniques  and  processes  were  used  in  test 
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specimen  fabrication  .  All  ribs  and  stringers  were  formed  from  flat  sheet  material  and  heat 
treated  as  necessary.  All  heat  treatment  was  checked  against  process  standards.  Quality 
control  procedures  were  used  at  all  times.  The  material  of  all  parts  for  the  test  specimens 
described  here  was  7075-T6. 

Each  specimen  was  assembled  using  the  test  frame  as  a  jig  to  insure  a  minimum  of  specimen 
prestress  in  the  test  configuration.  The  specimen  fabrication  sequence  was  such  that  the 
acoustic  fatigue  testing  was  continued  without  interruption . 

1 .  Stiffened  Panel  Specimens 

The  configuration  of  the  stiffened  panel  specimens  is  illustrated  in  Figure  15.  The  basic 
details  of  the  structure  are  given  in  Table  V.  The  stringer  dash  number  given  in  Table  V 
refers  to  the  stringer  detail  dimensions  given  in  Table  VI.  For  all  specimens,  the  stringers 
were  parallel  to  the  length  of  the  specimen  .  The  channel  section  rib  details  are  also 
presented  in  Table  VI.  Figure  16  illustrates  the  overall  specimen  test  configuration  and 
the  method  used  for  mounting  the  specimens  in  the  test  fixture.  Both  the  skin  side  and  the 
rib-stringer  side  are  illustrated.  The  specimen  designation  follows  the  numerical  sequence 
established  in  Reference  8. 


TABLE  V 

STIFFENED  PANEL  DESIGN  DETAILS 


Designation 

Qty . 

a 

b 

STR-31 

2 

6.0 

12.0 

STR-32 

2 

6.0 

18.0 

STR-33 

2 

6.0 

12.0 

STR-34 

2 

6.0 

18.0 

STR-35 

2 

6.0 

12.0 

STR-36 

2 

6.0 

18.0 

STR-37 

2 

9.0 

18.0 

STR-38 

2 

9.0 

18.0 

STR-39 

2 

6.0 

18.0 

STR-40 

2 

9.0 

18.0 

STR-41 

2 

6.0 

12.0 

STR-42 

2 

6.0 

12.0 

NOTE:  All  dimensions  in  inches. 
*See  Table  VI . 


t 

s 

Rib 

Stringer 

Stringer 

shape 

0.040 

-1* 

-1* 

zee 

0.040 

-1 

-1 

zee 

0.040 

-1 

-2 

channel 

0.040 

-1 

-2 

channel 

0.040 

-1 

-3 

hat 

0.040 

-l 

-3 

hat 

0.040 

-2 

-4 

zee 

0.040 

-2 

-5 

channel 

0.050 

-2 

-4 

zee 

0.032 

-2 

-1 

zee 

0.032 

-2 

-6 

zee 

0.032 

-2 

-7 

channel 
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CLIP  DETAIL:  ZEE  SECTION  STRINGER 


CLIP  DETAIL:  CHANNEL  SECTION  STRINGER 


CLIP  DETAIL:  HAT  SECTION  STRINGER 
FIGURE  15  STIFFENED  PANEL  SPECIMEN  CONFIGURATION  (CONCLUDED) 
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TABLE  VI 

STIFFENED  PANEL  STRINGER  DESIGN  DETAILS 
t  b,  b 


Stringer 

Shope 

r 

-1 

zee 

0.032 

-2 

channel 

0.032 

-3 

hat 

0.032 

-4 

zee 

0.040 

-5 

channel 

0.040 

-6 

zee 

0.025 

-7 

channel 

0.025 

t 

Rib 

Shape 

r 

-1 

channel 

0.080 

-2 

channel 

0.050 

NOTE:  All  dimensions  in  inches 


1.25 

0.75 

- 

0.09 

1.25 

0.75 

- 

0.09 

1 .00 

0.50 

0.75 

0.09 

1 .25 

0.75 

- 

0.09 

1.25 

0.75 

- 

0.09 

1.25 

0.75 

- 

0.09 

1.25 

0.75 

- 

0.09 

h 

r 

2.90 

0.75 

0.22 

2.90 

0.75 

0.25 

-*l  b1  U- 


BRACE 


2.  Nine  Cell  Box  Structure 


The  configuration  of  the  box  specimen  structure  is  illustrated  in  Figure  17.  The  details  of 
the  structure  are  given  in  Table  VII .  For  purpose  of  discussion,  the  internal  structure  shall 
be  classified  as  either  a  frame  or  a  rib.  For  all  specimen  designs,  the  frame  was  continuous 
across  the  width  of  the  specimen .  The  ribs  were  divided  into  a  center  rib  (test  rib)  and  an 
end  rib.  The  center  rib  and  the  end  rib  were  attached  to  the  frames  with  angle  clips  along 
the  height  of  the  specimen  (dimension  h) .  An  angle  clip  was  attached  to  the  end  of  each 
frame  and  rib  for  attaching  the  specimen  to  the  test  fixture.  The  frames  and  the  ribs  were 
oriented  so  that  the  flanges  faced  the  edge  of  the  specimen.  The  test  fixture,  also  used  as 
an  assembly  jig,  was  designed  to  allow  for  variation  of  frame  and  rib  spacing  and  specimen 
height.  The  cover  sheets  and  ribs  were  firmly  bolted  to  the  test  frame  along  all  edges. 

This  fixture  design  allowed  visual  inspection  of  all  structure  surfaces  in  the  test  configura¬ 
tion  .  Figure  18  illustrates  the  details  of  the  test  frame  as  viewed  from  the  rear  with  test 
specimen  BX6B  installed. 


C.  General  Test  Procedure 


The  various  tests  were  conducted  using  the  following  procedure: 

1 .  Modal  Frequency  Studies 

1 .1  Stiffened  Panel  Specimens  -  The  pair  of  stiffened  panel  specimens,  mounted  in  the 
test  frame  as  illustrated  in  Figure  16,  were  placed  over  an  enclosure  containing  two 
electromechanical  speakers  as  illustrated  in  Figure  19  with  the  stiffeners  exposed  for 
observation .  Cork  particles  were  sprinkled  on  the  test  specimen  and  a  frequency  sweep 
conducted .  The  Chladni  patterns  formed  by  the  cork  particles  were  observed  and  sketched . 
For  the  most  predominant  modes,  accelerometer  measurements  were  conducted  to  determine 
panel  and  stiffener  motion  . 

1.2  Box  Specimens  -  The  modal  frequency  studies  for  the  box  specimens  were  conducted 
in  a  different  manner  than  that  des~;ibed  for  the  stiffened  panel  specimens.  With  the 
specimen  mounted  in  the  test  frame  as  illustrated  in  Figure  18,  two  speaker  enclosures 
were  placed  on  either  side  of  the  specimen  as  illustrated  in  Figure  20.  Accelerometers 
were  placed  on  the  center  of  each  panel  between  the  frames  and  the  test  rib.  A  frequency 
sweep  was  conducted  and  the  accelerometer  output  recorded  on  an  x-y  plot.  Three 
speaker  phase  conditions  were  used  in  this  study.  For  purposes  of  discussion,  the  side  of 
the  structure  to  be  exposed  to  random  acoustic  excitation  will  be  called  the  front.  The 
speaker  excitation  conditions  used  in  the  modal  frequency  studies  were:  one  sided  excita¬ 
tion  speakers  on  the  front  side  only;  two  sided  excitation  speakers  on  both  sides,  in  phase 
and  out-of-phase.  Cork  particles  were  sprinkled  on  the  test  ribs,  and  the  Chladni  patterns 
formed  for  each  mode  were  sketched.  Using  accelerometers,  the  phase  relations  between 
the  panel  bays  and  the  ribs  were  determined.  The  speaker  excitation  conditions  ensured 
that  all  possible  modes  were  excited.  To  determine  the  panel  modes  on  the  front  surface 
of  the  structure,  the  specimen/test  frame  were  placed  with  the  front  face  up  and  a  speaker 
enclosure  placed  at  an  angle  to  the  specimen  to  allow  visual  observation  of  the  specimen 
surface.  Cork  particles  were  sprinkled  over  the  surface  of  the  structure  and  Chladni 
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-Hh-  0.32  TYP 


UPPER  (TEST)  SURFACE  SHOWN 
PLAN  VIEW  SPECIMEN 


NOTE:  LOWER  (REAR)  SURFACE  RIVETS  STSAD004A05-(  ) 
UNIVERSAL  HEAD,  E.D.  =0.32 

FIGURE  17  BOX  SPECIMEN  CONFIGURATION 


FRAME  DETAIL 


I 

rr 

h 

ill 


END  RIB  DETAIL 


DETAIL  D 


h  • 


CENTER  RIB  DETAIL 


FIGURE  17  BOX  SPECIMEN  CONFIGURATION  (CONCLUDED) 
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TABLE  VII 


BOX  SPECIMEN  DESIGN  DETAILS 


Internal  Structure 


Cover  Sheet 

rib 

frame 

Designation 

Qty . 

a 

b 

t 

s 

h 

t 

r 

r 

!f  i 

BX-1 

2 

10.0 

20.0 

0.063 

10.0 

0.063 

0.16 

0.080  0.22 

BX-2 

2 

8.0 

24.0 

0.050 

10.0 

0.063 

0.16 

0.080  0.22 

BX-3 

2 

7.0 

28.0 

0.040 

10.0 

0.050 

0.12 

0.080  0.22 

BX-4 

2 

10.0 

20.0 

0.063 

10.0 

0.032 

0.16 

0.040  0.22 

BX-5 

2 

8.0 

24.0 

0.040 

10.0 

0.020 

0.12 

0.040  0.22 

BX-6 

2 

8.0 

24.0 

0.040 

7.0 

0.020 

0.12 

0.040  0.22 

NOTE:  All  dimensions  in  inches. 
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FIGURE  19  STIFFENED  PANEL  SPECIMEN  MOUNTED  FOR 
MODAL  FREQUENCY  STUDY 


(b)  SPECIMEN  CONFIGURATION  FOR  OBSERVING  FRONT 
COVER  SHEET  MODES 


FIGURE  20  BOX  SPECIMEN  MOUNTED  FOR  MODAL  FREQUENCY  STUDY 
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patterns  for  each  predominant  mode  were  sketched.  With  the  response  plots,  phase  infor 
motion,  and  Chladni  patterns,  it  was  possible  to  define  each  responding  mode. 

2.  Test  Specimen  Instrumentation 


2.1  Stiffened  Panel  Specimens  -  After  a  thorough  study  of  the  test  specimen  response 
characteristics,  strain  gage  locations  were  selected.  Figure  21  shows  typical  uniaxial 
strain  gage  locations  for  the  stiffened  panel  specimens. 

2.2  Box  Specimens  -  Uniaxial  strain  gage  installation  on  the  box  specimens  varied; 
however,  typical  installations  for  the  front  surface  and  ribs  is  indicated  in  Figure  22.  For 
the  speaker  excitation  conditions  described  in  Section  1  .2,  frequency  sweeps  for  selected 
strain  gages  were  conducted  with  the  strain  gage  response  recorded  on  an  x-y  plot. 

3.  Damping  Studies 

Frequency  sweeps  using  the  electromechanical  speakers  previously  described  were  con¬ 
ducted  with  the  strain  gage  response  displayed  on  an  oscilloscope.  At  frequencies  with 
predominant  strain  gage  response,  the  excitation  was  suddenly  chopped,  and  the  decaying 
strain  signal  on  the  oscilloscope  photographed.  The  decaying  strain  signals  so  obtained 
were  used  to  determine  the  logarithmic  decrement,  6,  or  the  damping  ratio,  £,  Reference 
14.  For  the  box  specimens,  only  excitation  on  the  front  side  was  used  for  the  damping 
studies. 

4.  High  Intensity  Sinusoidal  Frequency  Sweeps 

The  test  specimens,  instrumented  as  described  in  Section  2,  were  placed  on  the  progressive 
wave  test  section  of  the  high  intensity  acoustic  fatigue  facility  at  the  Aerospace  Sciences 
Research  Laboratory.  Typical  installation  details  are  described  in  Reference  8. 

Grazing  incidence  sinusoidal  frequency  sweeps  in  the  range  of  50-1000  Hz  were  conducted 
at  the  selected  specimen  test  level .  Strains  for  all  gages  were  recorded  on  x-y  plots.  The 
strain  response  was  studied  and  used  to  shape  the  broad-band  acoustical  excitation  for  the 
fatigue  tests. 

5.  Broad-Bond  Acoustical  Noise  Test  Spectra 

The  test  specimen/frame  assembly  was  removed  from  the  progressive  wave  test  section 
and  replaced  by  a  one-inch  thick  plywood  panel .  The  test  sound  pressure  spectrum  was 
then  shaped  out  of  the  presence  of  the  vibrating  test  specimen.  The  width  of  the  exci¬ 
tation  spectrum  was  at  least  three  times  the  width  of  the  strain  response  of  the  panel 
obtained  during  the  high  intensity  sinusoidal  sweeps.  Spectrum  shaping  was  required  to 
concentrate  the  acoustical  energy  in  the  desired  frequency  range.  The  plywood  panel 
was  removed  when  the  desired  spectrum  had  been  obtained. 

6.  Acoustic  Fatigue  Tests 

The  test  spec imen/Frame  assembly  was  re-installed  in  the  progressive  wave  test  section 
ready  for  fatigue  testing  to  take  place.  Test  specimen  inspection  was  made  in  the  following 
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SECTION 
A -A 


SECTION 

B-B 


SG9 


SGI  5 


SG8 


SG13— 

(b)  STRAIN  GAGE  LOCATIONS  FOR  TEST  RIB 


FIGURE  22  TYPICAL  STRAIN  GAGE  LOCATIONS  FOR  BOX  SPECIMENS 


sequence:  every  five  minutes  for  the  first  hour;  every  fifteen  minutes  for  the  next  two 
hours;  every  half  hour  for  the  next  seven  hours;  and  every  hour  until  test  termination. 

When  a  fatigue  failure  was  first  observed,  the  time,  location,  and  a  brief  description 
was  recorded  in  a  laboratory  record  book.  The  time  of  failure  and  crack  length  were 
marked  on  the  specimen.  Cracks  were  allowed  to  propagate  until  complete  modes  of 
failure  had  been  observed.  A  complete  description  of  the  data  recording  and  analysis 
systems  used  for  mode  studies  and  fatigue  tests  are  presented  in  Appendix  II. 


D.  Test  Results 


Results  of  the  stiffened  panel  and  box  specimen  experiments  are  presented  in  the  follow¬ 
ing  sections. 


1.  Stiffened  Panel  Test  Results 


Each  stiffened  panel  design  consisted  of  two  specimens.  Following  the  notation  used  by 
Ballentine,  the  specimen  nomenclature  is  given  by  the  panel  design  with  the  suffix  A 
and  B  to  denote  the  two  specimens.  For  all  fatigue  tests,  the  stiffened  panel  specimens 
A  and  B  were  tested  simultaneously. 


1.1  Modal  Frequency  Studies  -  As  discussed  in  Reference  5,  the  types  of  modes  to  be 
expected  for  nine  bay  panels  would  consist  of  various  phase  relations  between  adjacent 
bays  for  each  mode  number  (m,n)  of  a  panel .  The  Chladni  patterns  observed  during  the 
mode  studies  described  in  Section  0.1.1  are  presented  in  Appendix  III.  The  similarity 
of  the  predominant  modes  for  each  panel  design  is  apparent.  The  predominant  mode  for 
stiffener  response  was  a  fundamental  mode  for  each  panel  bay  between  ribs  with  adjacent 
bays  across  the  stringers  out  of  phase.  The  predominant  mode  for  skin  response  was  adja¬ 
cent  bays  across  ribs  out  of  phase  (reversed  bending).  Since  the  ribs  were  much  stiffer 
than  the  stringers,  the  strain  response  on  the  skin  was  generally  maximum  at  the  center 
of  the  short  side  of  the  center  bay.  This  was  to  be  expected  due  to  the  torsional ly  flex¬ 
ible  stringers  used  in  these  experiments.  A  summary  of  the  modal  frequencies  is  given 

in  Table  VIII. 

1.2  Damping  Ratios  -  The  decaying  strain  signals  as  described  in  Section  III.C.3  were 
used  to  obtain  average  values  for  the  damping  ratio.  An  arithmetic  average  of  the  damp¬ 
ing  ratio  for  each  design  was  obtained.  Table  IX  is  a  listing  of  these  average  damping 
ratios  as  determined  by  the  logarithmic  decrement  method.  There  were  no  significant 
differences  between  damping  for  the  skin  or  the  internal  structure.  Figure  23  presents  a 
plot  of  damping  ratio  versus  frequency  for  all  responding  strain  gages  and  all  specimens. 

1 .3  Hiah  Intensity  Frequency  Sweeps  -  Figures  24  and  25  ore  plots  of  strain  response  to 
hign  intensity  sinusoidal  excitation  at  the  indicated  sound  pressure  level.  The  frequency 
sweep  plots  are  typical  for  the  strain  responses  for  all  stiffened  panel  specimens  at  skin 
and  stringer  locations. 
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TABLE  VIII 

MODAL  FREQUENCIES  FOR  STIFFENED  PANELS* 
(Low  Intensity  Sinusoidal  Excitation) 


ecimen 

STR-31A 
STR-31A 
STR-31B 
STR-31B 
STR-32A 
STR-32B 
STR-33A 
STR-33B 
STR-34A 
STR-34A 
STR-34B 
STR-34B 
STR-35A 
STR-35B 
STR-36A 
STR-36A 
STR-36B 
STR-36B 
STR-37A 
STR-37B 
STR-38A 
STR-38B 
STR-39A 
STR-39B 
STR-40A 
STR-40B 
STR-41A 
STR-41B 
STR-42A 
STR-42B 

*Only  predominant  frequencies 


Frequency,  Kx 
160 

209,264 

165 

213,262 

166 
168 
194 
186 
159 
183 
158 
170 
205 
195 
166 
212 
180 
217 
75 
87 
96 
81 
177 
169 
52,92 
71 

123,164 
128,185 
180 
199 

listed  here,  see  Appendix  III  for 


Stringer  Motion 
Torsion 

Bending-torsion 

Torsion 

Bending-torsion 

Torsion 

Torsion 

Torsion 

Torsion 

Torsion 

Bending-torsion 

Torsion 

Bending-torsion 

Torsion 

Torsion 

Torsion 

Bending-torsion 

Torsion 

Bending-torsion 

Torsion 

Torsion 

Torsion 

Torsion 

Torsion 

Torsion 

Torsion 

Torsion 

Torsion 

Torsion 

Torsion 

Torsion 

Chladni  patterns. 
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TABLE  IX 


average  damping  ratios  for  stiffened  panels 

a - Damping  Ratio, _% 

Specimen 

— - 

2.3 

STR-31A 

1.8 

STR-31B 

2.8 

STR-32A 

2.8 

STR-32B 

2.1 

STR-33A 

2.0 

STR-33B 

1.9 

STR-34A 

1.3 

STR-34B 

2.1 

STR-35A 

2.5 

STR-35B 

1.4 

STR-36A 

1.8 

STR-36B 

1.4 

STR-37A 

1.5 

STR-37B 

1.1 

STR-38A 

2.1 

STR-38B 

1.0 

STR-39A 

1.3 

STR-39B 

1.9 

STR-40A 

1.5 

STR-40B 

1.8 

STR-41A 

1.0 

STR-41B 

1.5 

STR-42A 

1.0 

STR-42B 
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1 .4  Brood-Bond  Acoustical  Noise  Test  Spectra  -  The  acoustical  noise  test  spectra  were 
shaped  such  that  the  bandwidth  of  the  excitation  was  at  least  three  times  the  bandwidth  of 
the  test  specimen  sinusoidal  response.  If  more  than  one  mode  was  deemed  significant,  the 
test  spectrum  was  shaped  so  as  to  include  all  responding  modes.  Figures  26  and  27 are 
typical  broadband  test  spectra  for  the  stiffened  panel  designs.  All  narrow  band  analyses 
for  the  broadband  acoustical  noise  were  accomplished  using  a  10Hz  bandwidth  filter.  The 
test  spectrum  level  was  determined  using  the  simple  relation  Spectrum  Level  =  Sound 
Pressure  Level  -  10  log(&f),  dB,  where  Af  is  the  filter  bandwidth.  For  &f  =  10, 

Spectrum  Level  =  Sound  Pressure  Level  -  10  dB.  The  test  spectrum  level  is  indicated  in 
each  figure. 

1 .5  Fatigue  Tests  -  Twenty-four  stiffened  panel  specimens,  twelve  designs  with  two 
replicates  each,  were  exposed  to  broadband  random  acoustical  excitation  until  one  or 
more  cracks  developed  in  the  skin  or  stringer.  Table  X  is  a  summary  of  the  test  results. 

Magnetic  tape  recordings  were  made  of  the  signals  from  the  strain  gages  positioned  on 
the  skin  and  the  stringers.  These  strain  data  were  analyzed  using  a  narrow  band  filter 
(nominal  10  Hz)  to  determine  the  characteristics  of  the  strain  response  to  the  acoustical 
excitation.  Figures  28  through  35  are  the  narrow  band  analyses  of  typical  single  mode 
and  multi -mode  response. 

In  addition  to  the  narrow  band  analysis,  the  recorded  strain  data  was  analyzed  to  deter¬ 
mine  the  probability  density  of  instantaneous  strain.  Figures  36  through  43  are  typical 
results  of  these  statistical  analyses.  For  reference,  a  Gaussian  distribution  is  plotted  in 
each  figure. 

Figure  44  shows  the  relationship  between  overall  nominal  rms  stress  and  the  number  of 
cycles  to  failure  for  the  failure  data.  The  number  of  cycles  to  failure  was  established 
from  the  number  of  zero  crossings  with  negative  slope  (obtained  as  described  in  Appendix 
II)  for  each  recorded  strain  signal  and  the  time  failure  was  observed.  The  regression  line 
and  95%  confidence  limits  given  in  Figure  44  were  computed  as  described  in  Appendix  II, 
Reference  8.  The  highest  correlation  coefficient  for  these  data  was  obtained  using  failure 
data  for  strain  gages  located  on  the  stiffener  flange  at  the  attach  clip  to  the  frame  with 
the  strain  gage  axis  aligned  with  the  stringer  length.  The  regression  line  and  confidence 
limits  illustrated  <n  Figure  44  are  established  for  these  data  points  only  although  all  failure 
data  is  plotted . 

2.  Box  Specimen  Test  Results 

Due  to  the  physical  size  of  the  box  structure,  these  specimens  were  fatigue  tested  indi¬ 
vidually.  The  specimen  nomenclature  is  given  by  the  specimen  design  with  the  suffix  A 
and  B  to  denote  the  duplicate  specimens. 

2.1  Modal  Frequency  Studies  -  Mode  shapes  for  each  box  specimen  structure  were 
determined  as  described  in  Section  III  .C.l .  The  experimentally  observed  modal  frequencies 
are  given  in  Table  XI  for  each  specimen .  The  notation  for  a  specific  mode  is  given  by  the 
mode  numbers  (m,n,q)  as  described  in  Section  II.  Figures  45  through  53  illustrate  the 
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FIGURE  26  BROAD  BAND  TEST  SPECTRUM:  STR  37  A  &  B 
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FIGURE  27  BROAD  BAND  TEST  SPECTRUM:  STR  39A  &  B 
(MULTIMODE  RESPONSE) 
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TABLE  X 

SUMMARY  OF  STIFFENED  PANEL  FATIGUE  TESTS 


Specimen 

Test 

Spectrum 
Level,  dB 

Resp~  •  :e 
Frequency 
Hz 

N  x  10“6 
c 

a 

ksi 

Remarks 

STR-31A 

124 

160 

2.2 

5.2  (FT)* 

B 

124 

170 

- 

- 

No  failure  at  gage 

STR-32A 

124 

180 

5.5 

2.0  (FT) 

A 

124 

180 

6.8 

3.2  (FT) 

A 

124 

180 

5.2 

4.2  (FA) 

A 

124 

180 

5.1 

2.0  (WA) 

B 

124 

170 

- 

- 

No  failure  at  gage 

STR-33A 

124 

160,350 

8.1 

2.2  (FA) 

B 

124 

140,330 

6.8 

2.3  (FA) 

B 

124 

140,330 

7.8 

5.6  (FT) 

STR-34A 

124 

160,350 

5.2 

3.4  (WA) 

A 

124 

140,330 

8.9 

6.8  (FT) 

B 

124 

- 

- 

No  failure  at  gage 

STR-35A 

124 

240 

- 

- 

No  failure  at  gage 

B 

124 

220 

16.1 

1 .4  (FA) 

B 

124 

220 

16.7 

0.8  (FA) 

STR-36A 

124 

160 

- 

- 

No  failure  at  gage 

B 

124 

185 

8.3 

2.8  (FA) 

B 

124 

185 

7.0 

3.6  (FA) 

STR-37A 

132 

145 

4.2 

1 .6  (WA) 

B 

132 

140 

1.4 

1 .2  (WA) 

B 

132 

140 

2.2 

2.3  (WA) 

STR-38A 

128 

175 

2.2 

5.2  (WA) 

A 

128 

175 

0.7 

3.5  (FA) 

A 

128 

175 

1.03 

4.9  (FA) 

B 

128 

140 

2.5 

4.2  (FA) 
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TABLE  X  (Continued) 


Specimen 

Test 

Spectrum 
Level,  dB 

Response 

Frequency 

Hz 

N  x  10"6 

c 

0 

ksi 

Remarks 

STR-39A 

132 

200,370 

1.4 

4.4  (WA) 

A 

132 

200,370 

0.7 

4.5  (WA) 

B 

132 

200,400 

1.0 

4.2  (WA) 

B 

132 

200,400 

1.8 

3.7  (WA) 

STR-4QA 

132 

140,300 

1.22 

16.0  Skin  @  Frame 

A 

132 

140,300 

0.44 

8.7  (FA) 

A 

132 

140,300 

0.37 

6.7  (FA) 

A 

132 

140,300 

0.41 

12.9  (FA) 

A 

132 

140,300 

0.43 

8.8  (FA) 

B 

132 

140,300 

1.5 

3.6  (FA) 

B 

132 

140,300 

0.36 

4.3  (FA) 

B 

132 

140,300 

1.6 

4.3  (FA) 

B 

132 

140,300 

0.40 

4.6  (FA) 

STR-41A 

124 

220 

6.6 

2.1  (FA) 

A 

124 

220 

10.9 

2.6  (FA) 

B 

124 

220 

9.7 

2.3  (FA) 

STR-42A 

124 

200-230 

- 

2.1  (FA) 

No  failure  at  gage 

A 

124 

200-230 

- 

2.6  (FA) 

No  failure  at  gage 

A 

124 

200-230 

- 

3.4  (FA) 

No  failure  at  gage 

A 

124 

200-230 

- 

2.4  (FA) 

No  failure  at  gage 

B 

124 

200-230 

0.89 

3.2  (FA) 

B 

124 

200-230 

- 

3.8  (FA) 

No  failure  at  gage 

B 

124 

200-230 

- 

2.5  (FA) 

No  failure  at  gage 

B 

124 

200-230 

- 

2.1  (FA) 

No  failure  at  gage 

*(FT)  -  Flange  transverse  gage  at  stiffener  center 
(FA)  -  Flange  axial  gage  at  clip 
(WA)  -  Web  axial  gage  at  clip 
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SPECTRUM  LEVEL:  132  dl 
FILTER  BANDWIDTH:  10 
SCAN  RATE:  5.5  Hz/sec 
OVERALL  STRAIN:  e-  8 
LOCATION:  SEE  FIGUR 


narrow  BAND 


FREQUENCY,  Hz 

ANALYSIS:  STRAIN  GAGE  1 ,  STR-37A 


FIGURE  28 


RMS  STRAIN 
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SPECTRUM  LEVEL:  132  dB 
FILTER  BANDWIDTH:  10  Hz 
-f  SCAN  RATE:  5.5  Hz/sec 

OVERALL  STRAIN:  €  =  240 |j in/in 
LOCATION:  SEE  FIGURES  24  &  25 
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FIGURE  31  NARROW  BAND  ANALYSIS:  STRAIN  GAGE  6,  STR-37A 
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|  SPECTRUM  LEVEL:  132  dB 
FILTER  BANDWIDTH:  10  Hz 
SCAN  RATE:  5.5  Hz/sec 
OVERALL  STRAIN:  c=  200u  in/in  I 
LOCATION:  SEE  FIGURE  25 


FREQUENCY,  Hz 


FIGURE  32  NARROW  BAND  ANALYSIS:  STRAIN  GAGE  7,  STR-37A 
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FIGURE  33  NARROW  BAND  ANALYSIS:  STRAIN  GAGE  1,  STR-39A 
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RMS  STRAIN 


FREQUENCY,  Hi 

FIGURE  35  NARROW  BAND  ANALYSIS:  STRAIN  GAGE  6,  STR-39A 
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INSTANTANEOUS  TO  RMS  RATIO 


FIGURE  38  AMPLITUDE  DISTRIBUTION  OF  STRAIN:  STRAIN  GAGE  6,  STR-37A 
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INSTANTANEOUS  TO  RMS  RATIO 

FIGURE  39  AMPLITUDE  DISTRIBUTION  OF  STRAIN:  STRAIN  GAGE  9,  STK-37A 
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INSTAhttANtOUS  TO  RMS  RATIO 

FIGURE  40  AMPLITUDE  DISTRIBUTION  OF  STRAIN:  STRAIN  GAGE  1,  STR-39A 
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INSTANTANEOUS  TO  RMS  RATIO 

FIGURE  41  AMPLITUDE  DISTRIBUTION  OF  STRAIN:  STRAIN  GAGE  3,  STR-39B 


FIGURE  42  AMPLITUDE  DISTRIBUTION  OF  STRAIN:  STRAIN  GAGE  6,  STR-39A 
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INSTANTANEOUS  TO  RMS  RATIO 


FIGURE  43  AMPLITUDE  DISTRIBUTION  OF  STRAIN:  STRAIN  GAGE  9,  STR-39A 
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FIGURE  44  STIFFENED  PANEL  STRINGER  FATIGUE  CURVE 


TABLE  XI 


MODAL  FREQUENCIES  FOR  BOX  SPECIMENS 
(Low  Intensity  Sinusoidal  Excitation) 


Mode 


Specimen 

0,1,1) 

0,2,1) 

0,3,1) 

0,3, 

BX-1A 

94 

191 

222 

— 

BX-1B 

78 

144 

235 

— 

BX-2A 

90 

136 

193 

— 

BX-2B 

92 

135 

220 

— 

BX-3A 

82 

— 

137 

— 

BX-3B 

81 

106 

133 

— 

BX-4A 

139 

181 

228 

— 

BX-4B 

135 

— 

209 

— 

BX-5A 

78 

92 

164 

134 

BX-5B 

88 

91 

— 

136 

BX-6A 

62 

77 

— 

124 

BX-6B 

62 

— 

— 

116 
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ACCELERATION,  g’s 


FREQUENCY,  Hz 

FIGURE  50  CENTER  RIB  RESPONSE:  OUT  OF  PHASE  EXCITATION  BOTH  SIDES,  BX4B 
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ACCELERATION,  g's 


acceleration  response  of  the  front  center  panel  bay,  the  center  rib,  and  the  rear  center 
panel  bay  for  the  three  excitation  conditions  of  the  low  intensity  sine  sweeps.  The 
accelerometer  location  was  the  center  of  the  bay  and  rib.  These  results  are  typical  of  the 
low  intensity  sine  sweeps  conducted  during  the  modal  frequency  studies. 

2.2  Damping  Ratios  -  The  decaying  strain  signals  as  described  in  Section  III  .C  .3  were 
used  to  obtain  average  values  for  the  damping  ratio.  An  arithmetic  average  for  each 
design  was  obtained.  Table  XII  is  a  listing  of  these  average  damping  ratios  as  determined 
by  the  logarithmic  decrement  method.  A  plot  of  damping  ratio  versus  frequency  for  the 
box  specimens  is  given  in  Figure  23. 

2.3  High  Intensity  Frequency  Sweeps  -  Figures  54  through  57 are  plots  of  strain  response 
to  high  intensity  sinusoidal  excitation  at  the  indicated  sound  pressure  level.  The  strain 
gages  illustrated  are  for  specimen  BX4B.  Strain  gage  1  is  located  across  the  rivet  line  at 
the  middle  of  the  long  side  of  the  center  bay  on  the  side  exposed  to  acoustical  excitation. 
Strain  gage  7  is  located  across  the  rib  flange  directl>  under  gage  1 .  Strain  gage  8  is 
located  in  the  center  of  the  rib  with  the  gage  axis  in  the  direction  of  the  rib  height. 

Strain  gage  9  is  located  across  the  rib  flange  directly  opposite  to  strain  gage  7.  These 
figures  illustrate  the  discrete  frequency  strain  response  across  the  height  of  the  structure 
with  excitation  on  one  side. 

2.4  Broadband  Acoustical  Noise  Test  Spectra  -  The  acoustical  noise  test  spectra  were 
shaped  such  that  the  bandwidth  of  the  excitation  was  at  least  three  times  the  bandwidth  of 
the  test  specimen  sinusoidal  response.  If  more  than  one  mode  appeared  to  be  significant, 
the  test  spectrum  was  shaped  so  as  to  include  all  responding  modes.  Figures  58  through  60 
illustrate  typical  broadband  test  spectra  for  the  box  specimens.  Figure  60  illustrates  the 
test  spectrum  for  specimen  BX4B. 

2.5  Fatigue  Tests  -  Twelve  box  structure  specimens,  six  designs  with  two  replicates  each, 
were  exposed  to  broadband  random  acoustical  excitation  until  one  or  more  cracks  developed 
in  the  skin  or  ribs.  Table  XIII  is  a  summary  of  the  test  results. 

Magnetic  tape  recordings  were  made  of  the  signals  from  the  strain  gages  positioned  in  the 
skin  and  the  ribs.  These  strain  data  were  analyzed  using  a  narrow  band  filter  (nominal 
10  Hz)  to  determine  the  characteristics  of  the  strain  response  to  the  acoustical  excitation. 
Figures  61  through  64  are  narrow  band  analyses  of  the  strain  response  for  strain  gages  1 ,  7, 
8,  and  9  of  specimen  BX4B.  (See  Figures  54  through  57  for  the  high  intensity  sinusoidal 
response  of  these  strain  gages .) 

In  addition  to  the  narrow  band  analysis,  the  recorded  strain  data  were  analyzed  to  deter¬ 
mine  the  probability  density  of  instantaneous  strain.  Figures  65  through  68  are  plots  of 
amplitude  distribution  for  strain  gages  1 ,  7,  8,  and  9  of  specimen  BX4B.  These  results  are 
typical  of  the  statistical  analyses.  For  reference,  a  Gaussian  distribution  is  plotted  on 
each  figure. 

Figure  69  shows  the  relationship  between  overall  nominal  rms  stress  and  the  number  of 
cycles  to  failure  for  all  failure  data  .  The  number  of  cycles  to  failure  was  established  from 
the  number  of  zero  crossings  with  negative  slope  (obtained  as  described  in  Appendix  II)  for 


TABLE  XII 

AVERAGE  DAMPING  RATIOS  FOR  BOX  SPECIMENS 


Specimen 

Average  Damp! 

BX-1A 

1.7 

BX-1B 

2.5 

BX-2A 

2.2 

BX-2B 

2.0 

BX-3A 

2.4 

BX-3B 

2.2 

BX-4A 

1.3 

BX-4B 

1.4 

BX-5A 

2.2 

BX-5B 

3.0 

BX-6A 

1.5 

BX-6B 

2.0 

84 


FIGURE  54  HIGH  INTENSITY  SINE  SWEEP:  STRAIN  GAGE  1,  BX-4B 
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SPECTRUM  LEVEL,  dB  (Re:  0.0002  dyn/< 


TABLE  XIII 


SUMM 

Test 

Spectrum 

Level 

dB 

ARY  OF  BOX  SPECIMEN  FATIGUE  TESTS 

Response 

FT“’cy  N  x  !0'6 

Hz  c  ksi 

Mode 

Responding 

124 

140-150* 

6.2 

6.0  (SF)* 

(1/1,1) 

124 

140-150 

6.3 

6.8  (SF) 

0,1,1) 

124 

140-160 

11.5 

5.0  (SR) 

0,1,1) 

121 

140-180 

9.25 

5.2  (SF) 

(1,1,1)  &  (1,3,1) 

121 

100-160 

43.5 

6.6  (SR) 

(1,1,1)  &  (1,3,1) 

119 

130/160 

No  Failures 

At  Gages 

(1,1,1)  &  (1,3,1) 

119 

115-145 

No  Failures 

At  Gages 

(1,1,1)  40,3,1) 

124 

140-160 

4.3 

6.6  (UF) 

0,1,1) 

124 

140-160 

3.5 

7.0  (UF) 

0,1,1) 

124 

140-160 

4.2 

7.2  (UF) 

0,1,1) 

124 

140-160 

3.7 

5.7  (UF) 

0,1,0 

125 

150 

0.68 

13.9  (UF) 

0,1,1) 

125 

150 

0.94 

10.8  (UF) 

0,1,1) 

125 

150 

1.22 

10.3  (UF) 

(1,1,1) 

125 

150 

1.07 

11.3  (UF) 

0,1,1) 

124 

33/175 

0.22 

22.2  (UF) 

(1,1,1)40,3,2) 

124 

33/175 

0.25 

21.1  (UF) 

(1,1,1)40,3,2) 

124 

33/175 

0.25 

18.6  (UF) 

(1,1,1)40,3,2) 

124 

33/175 

0.23 

22.6  (UF) 

(1,1,1)40,3,2) 

125 

W/135 

0.25 

24.8  (UF) 

(1,1,1)40,3,2) 

125 

to/135 

0.77 

15.5  (LF) 

(1,1,1)40,3,2) 

125 

to/135 

0.28 

21 .6  (UF) 

(1,1,1)40,3,2) 

125 

to/135 

0.74 

16.5  (LF) 

(1,1,1)  4(1,3,2) 

125 

to/135 

0.27 

17.6  (UF) 

(1,1,1)40,3,2) 

125 

to/135 

0.29 

19.6  (UF) 

(1,1,1)40,3,2) 

I 


TABLE  XIII  (Continued) 


Test 


Specimen 

Spectrum 

Level 

dB 

Response 

Frequency 

Hz 

N  x  10"6 
c 

a 

ksi 

Mode 

Responding 

BX-5B 

125 

90/135* 

0.34 

19.6  (UF)* 

(1,1,1)  &  (1,3,2) 

BX-5B 

125 

90/135 

0.32 

18.6  (LF) 

(1,1,1)  &  (1,3,2) 

BX-6A 

124 

90/135 

0.47 

15.5  (UF) 

(1,1,1)  &  (1,3,1) 

BX-6A 

124 

90/135 

0.36 

16.5  (UF) 

(1,1,1)  &  (1,3,1) 

BX-6A 

124 

90/135 

0.55 

13.9  (UF) 

(1,1,1)  &  (1,3,1) 

BX-6A 

124 

90/135 

1.78 

7.4  (LF) 

(1,1,1)  &  (1,3,1) 

BX-6A 

124 

90/135 

2.50 

9.8  (LF) 

(1,1,1)  &  (1,3,1) 

BX-6A 

124 

90/135 

2.30 

9.1  (LF) 

(1,1,1)  &  (1,3,1) 

BX-6A 

124 

90/135 

0.84 

11.3  (UF) 

(1,1,1)  &  (1,3,1) 

BX-6A 

124 

90/135 

0.95 

12.4  (UF) 

(1,1,1)  &  (1,3,1) 

BX-6A 

124 

90/135 

1.14 

11.9  (UF) 

(1,1,1)  &  (1,3,1) 

BX-6A 

124 

90/135 

2.24 

8.5  (LF) 

0,1,1)  Ml, 3,1) 

BX-6A 

124 

90/135 

2.10 

7.0  (LF) 

(1,1,1)  4(1, 3,1) 

BX-6B 

124 

95/125 

0.16 

18.5  (UF) 

(1,1,1)40,3,1) 

BX-6B 

124 

95/125 

0.34 

13.9  (UF) 

(1,1,1)  &(1,3,1) 

BX-6B 

124 

95/125 

0.31 

16.0  (UF) 

(1,1,1)  &  (1,3,1) 

BX-6B 

124 

95/125 

0.79 

16.5  (LF) 

(1,1,1)  &(1,3,1) 

BX-6B 

124 

95/125 

1.02 

12.9  (LF) 

(1,1,1)40,3,1) 

BX-6B 

124 

95/125 

0.89 

13.4  (LF) 

(1,1,1)  4(1,3,1) 

BX-6B 

124 

95/125 

0.50 

18.0  (UF) 

(1,1,1)40,3,1) 

BX-6B 

124 

95/125 

0.58 

16.0  (UF) 

(1,1,1)  4(1,3,1) 

BX-6B 

124 

95/125 

0.84 

14.9  (UF) 

(1,1,1)40,3,1) 

BX-6B 

124 

95/125 

0.50 

20.6  (UF) 

(1,1,1)  4(1,3,1) 

BX-6B 

124 

95/125 

0.52 

15.0  (LF) 

(1,1,1)40,3,1) 

BX-6B 

124 

95/125 

1.07 

10.3  (LF) 

(1,1,1)  4(1,3,1) 

"Notes:  140-150  denotes  a  broad  response  SF  -strain  gage  on  skin  across  frame 
in  frequency  range .  SR  -  strain  goge  on  skin  across  rib 

83/175  denotes  a  multimode  UF- strain  gage  on  rib  flange,  excited  surface 

response  at  indicted  LF  -strain  gage  on  rib  flange,  unexcited  surface 

frequenc  ies. 
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RMS  STRAIN,  in/in/^Hz  x  10 
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CM 


zt  SPECTRUM  LEVEL:  125  dB 
FILTER  BANDWIDTH:  10  Hz 
SCAN  RATE:  5.5  Hz/sec. 

— f-  OVERALL  STRAIN:?  =  840  u  in/ 
.  LOCATION:  See  Figure  22 
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FREQUENCY,  Hz 


FIGURE  61  NARROW  BAND  STRAIN  ANALYSIS:  STRAIN  GAGE  1 ,  BX-4B 


RMS  STRAIN 


FREQUENCY,  Hz 

FIGURE  64  NARROW  BAND  STRAIN  ANALYSIS:  STRAIN  GAGE  9,  BX-4B 


INSTANTANEOUS  TO  RMS  RATIO 


INSTANTANEOUS  TO  RMS  RATIO 


FIGURE  66  AMPLITUDE  DISTRIBUTION  OF  STRAIN:  STRAIN  GAGE  7,  BX-4B 
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INSTANTANEOUS  TO  RMS  RATIO 

FIGURE  67  AMPLITUDE  DISTRIBUTION  OF  STRAIN:  STRAIN  GAGE  9,  BX-4B 


INSTANTANEOUS  TO  RMS  RATIO 


FIGURE  68  AMPLITUDE  DISTRIBUTION  OF  STRAIN:  STRAIN  GAGE  8,  BX-4B 
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each  recorded  strain  signal  and  the  time  failure  was  observed.  The  regression  line  and 
95%  confidence  limits  given  in  Figure  69  were  computed  as  described  in  Appendix  II, 
Reference  8.  The  highest  correlation  coefficient  for  these  data  was  obtained  by  consid¬ 
ering  only  failure  data  for  strain  gages  located  on  the  rib  flange  across  the  rivet  line. 
The  regression  line  and  confidence  limits  illustrated  in  Figure  69  are  established  for 
these  data  only,  although  all  failure  data  is  plotted. 


E.  Discussion  of  Test  Results 


1.  Stiffened  Panel  Results 


1.1  Modal  Frequency  Studies  -  For  nine  bay  stiffened  panels  with  the  stringers  much 
more  flexible  in  bending  and  torsion  than  the  ribs,  it  was  expected  that  the  center  panel 
row  between  ribs  {test  area)  would  be  the  most  responsive.  This  was,  generally,  the  case. 
The  predominant  mode  was  a  fundamental  in  each  panel  bay  with  adjacent  bays  out-of- 
phase  (stringer  torsion).  The  stringer  bending  mode  (adjacent  bays  in  phase)  was  observed 
as  indicated  in  Table  VIII;  however,  these  modes  were  not  so  predominant  as  the  stringer 
torsion  modes.  The  reason  for  similar  modes  at  different  frequencies  is  due  to  the  degree 
of  coupling  between  adjacent  bays  influenced  by  the  stringer  bending  -  torsion  charac¬ 
teristics.  Acceleration  measurements  confirmed  that  the  stringer  deformation  was  basic¬ 
ally  a  twist  for  all  cross-section  shapes  considered. 

A  comparison  of  the  frequencies  determined  in  the  mode  study  and  the  response  frequen¬ 
cies  determined  from  the  high  intensity  size  sweeps  (Table  X)  show  significant  increases 
in  modal  frequency  for  the  fatigue  test  configuration.  This  increase  in  frequency  is  at¬ 
tributed  to  temperature  decrease  (as  much  as  40°F)  between  the  room  in  which  the  mode 
studies  were  conducted  and  the  acoustic  fatigue  facility. 

1.2  Frequency  Sweeps  -  The  high  intensity  sine  sweeps  were  conducted  to  determine 
the  specimen  response  in  the  fatigue  test  configuration.  For  ail  specimens  the  frequency 
range  used  to  establish  the  excitation  spectrum  was  based  upon  maximum  strain  response 
for  strain  gages  located  on  the  stringers.  Strain  gages  on  the  skin  generally  had  a  sim¬ 
ilar  frequency  response  as  gages  on  the  stringer  (see  Figure  24).  However,  if  the  mode 
study  indicated  that  the  panel  bays  across  a  frame  exhibited  appreciable  out-of-phase 
response  (reversed  bending)  a  strain  gage  was  located  on  the  frame  rivet  line  at  the  cen¬ 
ter  of  the  span.  These  strain  gages  experienced  a  multi-mode  response.  High  strain 
levels  at  the  frame  rivet  line  (short  siae  of  the  center  panel  row)  were  experienced  for 
both  the  fundamental  mode  of  the  center  panel  row  ana  the  mode  for  reversed  bending 
of  the  skin  across  the  frame. 

1 .3  Fatigue  Tests  -  Two  types  of  skin  failure  were  commonly  experienced  during  the 
fatigue  testing  of  the  stiffened  panel  specimens:  the  formation  of  fatigue  cracks  along 
a  skin  rivet  line  and  fatigue  failure  of  the  fasteners  (popped  rivet  heads).  These  types 
of  failure  occurred  only  along  the  frame  rivet  line  for  tne  specimens  tested. 

The  failure  experienced  by  design  STR-31  influenced  ail  subsequent  panel  designs.  Fig¬ 
ure  70  illustrates  the  fatigue  failure  mode  of  specimen  STR-31  A.  For  this  design,  the 
outer  bays  of  the  center  panel  row  were  of  greater  span  than  the  center  bay.  The  first 
observea  failure  occurred  at  the  center  of  me  rivet  line  on  the  short  side  (frame).  The 
ends  of  the  crack  were  drilled  and  a  cold  bonded  metal  patch  was  applied  to  the  skin 
side.  During  subsequent  testing,  the  cracks  continued  to  propagate  and  the  patch  de¬ 
laminated.  The  crack  propagated  very  quickly  along  the  edge  of  the  panel  and  into  the 
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center  bay  as  illustrated.  When  the  crack  passed  through  the  stringer,  the  flange  failed 
as  indicated  in  Figure  70.  This  failure  was  not  considered  to  be  va’id  since  the  failure 
did  not  originate  in  the  center  bay  test  area.  To  avoid  this  type  of  failure,  it  was  de¬ 
cided  that  the  center  bay  of  the  structure  should  be  equal  to  or  greater  in  area  than  the 
other  bays.  Since  the  total  specimen  size  was  determined  by  the  test  fixture,  the  cen¬ 
ter  bay  aspect  ratio  was  a  primary  variable,  and  the  nine  bay  configuration  was  desired, 
an  angle  brace  was  installed  along  the  length  of  the  outer  bays  of  the  center  panel  row. 
The  brace  comprised  of  an  aluminum  angle  attached  to  the  cage  pan  and  riveted  to  the 
skin.  This  design  change  decreased  the  span  of  the  outer  bays  so  that  the  center  panel 
bay  had  the  greatest  surface  area.  The  brace  is  illustrated  in  Figure  15.  It  was  in¬ 
stalled  on  all  designs  except  STR-37,  -38,  and  -40  (see  dimension  'a1  Table  V). 

Two  types  of  stringer  failures  were  experienced  for  the  specimens  tested  in  this  program: 
failure  at  the  clip  attachment  to  the  frame  and  failure  at  the  flange  rivet  line  along  the 
center  of  the  stringer  span.  All  designs  tested  experienced  failures  at  the  clip  attach¬ 
ment  to  the  frame.  Only  channel  section  stiffeners  experienced  fatigue  failures  at  the 
mid-span  of  the  flange  rivet  line.  Figure  71  illustrates  a  typical  clip  attachment  for 
zee  section  stringers:  two  rivets  along  the  web  height  attacking  the  stringer  to  the  clip. 
Channel  section  clip  attachments  were  similar.  The  clip  attachments  for  the  hat  section 
stiffeners  used  a  right  angle  clip  as  illustrated  in  Figure  15. 

Typical  stiffener  failures  at  the  clip  attachment  are  illustrated  in  Figure  72.  Figure  72a 
is  an  overall  view  of  the  frame-stringer  intersection  from  the  skin  side  (skin  removed) 
illustrating  the  failure  mode.  Typically,  a  crack  would  initiate  at  the  rivet  near  the 
upper  flange  and  propagate  down  the  bend  radius  and  across  the  upper  flange.  Figures 
72b  and  72c  show  closeup  views  of  such  failures  for  a  channel  section.  The  slight  bow 
in  the  upper  flange,  as  illustrated  in  Figure  72c,  was  apparently  caused  by  local  buck¬ 
ling  of  the  flange  after  fatigue  failure  was  complete.  For  this  specimen,  strain  gages 
were  installed  on  the  web  of  the  stringer  aligned  with  the  upper  clip  rivet.  The  strain  gage 
axis  is  mounted  parallel  to  the  stringer  axis  on  the  side  opposite  to  that  illustrated  (web 
axial  gage  (WA)  -  see  Figure  21  and  Table  X).  The  center  of  the  strain  gage  axis  was 
approximately  0.5"  in  from  the  point  at  which  the  crack  initiated.  Figure  72d  illus¬ 
trates  a  typical  stringer  failure  for  a  zee  section  for  which  the  crack  was  allowed  to 
propagate  down  the  stringer  web  and  flange  rivet  line.  The  strain  gage  installation  at 
this  failure  was  on  the  web  axis  between  the  rivet  heads  (overall  rms  stress  2.0  ksi).  A 
strain  gage  mounted  on  the  upper  flange  centered  on  the  clip  rivet  line  measured  an 
overall  rms  stress  of  4.0  ksi.  Both  of  these  strain  gages  were  within  0.25"  of  the  initial 
crack  location.  In  Figure  72d,  a  portion  of  the  flange  has  broken  off  (during  testing) 
allowing  the  complete  stringer  failure  to  be  viewed.  Figures  71  and  72  illustrate  the 
physical  impossibility  of  locating  a  strain  gage  at  the  initial  failure  point;  however,  as 
indicated  in  Figure  44,  the  failure  data  for  strain  gages  located  on  the  flange  and  aligned 
with  the  stringer  axis  yielded  the  most  significant  data  correlation.  The  measurement  of 
the  very  steep  strain  gradients  in  the  localized  area  of  the  stringer  clip  is  a  physically 
difficult  problem  that  can  only  be  resolved  by  additional  testing. 

The  stringer  flange  failures  at  mid-span  for  the  channel  section  are  illustrated  in  Fig¬ 
ure  73.  This  failure  occurred  independently  of  the  failures  illustrated  in  Figure  72a 
through  72c.  Figure  73a  is  an  overall  view  of  the  rib/frame  assembly  from  the  stiffener 
side  (skin  removed)  and  Figure  73b  is  a  closeup  view  of  the  failure.  In  Figure  73b  it  is 
interesting  to  note  that  cracks  initiated  at  the  edge  of  the  bucked  rivet  head  rather  than 
the  rivet  hole.  The  diagonal  failure  pattern  is  evidence  of  the  complex  strain  field  at 
this  location  on  the  stringer.  The  strain  gage  location  is  evident.  It  is  believed  that 
this  type  of  failure  was  due  to  stiffener  cross-section  distortion  in  the  plane  of  the  stif¬ 
fener  cross-soction .  This  type  of  failure  was  also  observed  on  the  flanges  of  the  box- 
specimen  structure. 
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FIGURE  71  TYPICAL  STRINGER  CLIP  ATTACHMENT 
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(b)  STRINGER  BENT 

I?J£P0SE  CRACK 
(OPPOSITE  END) 


(c)  CLOSEUP  VIEW  OF  FAILURE 


FIGURE  72  typical  stringer  failures  at  clip 
attachments 
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(a)  OVERALL  VIEW  OF  SPECIMEN 


(b)  CLOSEUP  VIEW  OF  FLANGE  FAILURE 
(INDEPENDENT  OF  FAILURES  AT  CLIPS) 


FIGUre  73  FLANGE  FAILURE  OF  CHANNEL  SECTION 
STRINGER 


Isolated  stringer/frame  clip  fatigue  failures  were  also  observed  during  testing;  however, 
these  failures  occurred  after  failure  of  the  stringer.  No  fatigue  failures  of  the  ribs  were 
recorded . 

2.  Box  Specimen  Structure 

2.1  Modal  Frequency  Studies  -  The  box  specimen  structure  exhibited  complex  dynamic 
behavior.  The  analytical  results  described  in  Section  II  were  an  aid  in  the  experimental 
investigation  in  that  a  uniform  mode  classification  was  developed.  As  indicated  in 
Table  XI,  only  a  few  modes  could  be  excited  and  identified  explicitly  as  coupled  skin- 
rib  modes.  The  degree  of  coupling  between  the  skin  and  rib  is  proportional  to  the  cube 
of  the  ratio  of  rib  to  skin  thickness  (see  Equations  39  and  45).  From  Table  VII,  it  can 
be  seen  that  designs  BX-1,  -2,  and  -3  represent  relatively  thick  ribs  whereas  designs 
BX-4,  -5,  and  -6  represent  relatively  thin  ribs  (as  compared  to  the  skin  thickness).  De¬ 
signs  BX-1,  -2,  and  -3  exhibited  a  higher  degree  of  coupling  between  the  skin  and  ribs 
than  did  designs  BX-4,  -5,  and  -6.  The  acceleration  response  of  the  cover  sheet  (front 
and  rear)  and  ribs  for  designs  BX-1,  BX-2,  and  BX-3  was  more  highly  coupled  than  that 
illustrated  in  Figures  45  through  53.  Although  the  theory  predicts  that  only  odd  number 
rib  modes  (across  the  structure  height)  would  be  excited  by  the  out-of-phase  condition 
(Eauation  46),  all  modes  were  excited  to  some  extent  for  all  speaker  conditions.  The 
only  significant  difference  in  response  was  observed  for  the  rear  cover  sheet  with  exci¬ 
tation  on  the  front  side.  As  would  be  expected,  the  rear  surface  response  was  signifi¬ 
cantly  less  (one  half  to  one  third)  than  that  measured  on  the  front  cover  sheet  for  one¬ 
sided  excitation  (see  Figure  51). 

2.2  Frequency  Sweeps  -  For  each  specimen,high  intensity  sine  sweeps  were  conducted 
to  determine  the  specimen  response  in  the  fatigue  test  condition.  For  all  specimens,  the 
frequency  range  used  to  establish  the  excitation  spectrum  was  based  upon  maximum  strain 
response  for  strain  gages  located  on  the  rib.  For  designs  BX-1,  -2,  and  -3  maximum  rib 
strain  response  occurred  at  the  center  of  the  rib  (strain  gage  8,  Figure  22).  For  designs 
BX-4,  -5,  and  -6  maximum  rib  strain  response  occurred  at  the  rib  flange  next  to  the 
front  surface  (exposed  to  excitation).  Strain  gages  located  across  the  frame  rivet  line 
exhibited  very  high  response  due  to  reversed  bending  of  the  skin  across  the  frame  (see 
Table  VII).  All  strain  gages  exhibited  predominant  response  in  the  frequency  intervals 
listed  in  Table  XIII.  Comparing  the  modal  frequencies  listed  in  Table  XI  with  the  re¬ 
sponse  frequencies  given  in  Table  XIII,  the  general  increase  in  modal  frequencies  indi¬ 
cated  is  attributed  to  the  temperature  difference  between  the  two  rooms  in  which  the 
tests  were  conducted  os  described  for  the  stiffened  panel  specimens. 

2.3  Fatigue  Tests  -  For  acoustic  fatigue  testing  the  box  specimen  structure  was  mounted 
in  the  test  section  side  wall  so  that  only  one  side  (front)  was  exposed  directly  to  the  ran¬ 
dom  noise.  For  all  specimens  skin  failures  were  observed  only  on  the  front  surface.  Two 
types  of  skin  failure  were  commonly  experienced  during  the  fatigue  testing.of  box  struc¬ 
ture:  the  formation  of  fatigue  cracks  along  a  skin  rivet  line  ana  fatigue  failure  of  the 
fasteners  (cracked  and/or  popped  rivet  heads).  The  skin  failures  occurred  only  for  the 
designs  BX-1,  -2,  and  -3.  Figure  74  illustrates  typical  skin  failures  for  specimens 
BX-1  A,  BX-1B,  BX-2A,  and  BX-3B.  The  location  and  extent  of  the  skin  failure  is  out¬ 
lined  for  each  specimen.  No  rib  failures  were  observed  for  designs  BX-1 ,  -2,  and  -3. 

Rib  fatigue  failures  were  observed  for  designs  BX-4,  BX-5,  and  BX-6.  No  skin  fatigue 
failures  for  these  designs  were  observed.  The  type  of  rib  failure  observed  is  illustrated 
in  Figure  75  for  specimen  BX-5B.  Figure  75a  illustrates  the  complete  specimen  with  the 
front  skin  removed.  The  fatigue  failure  of  the  ribs  was  complete  along,  the  entire  rib. 
length.  Figure  75b  is  a  closeup  view  of  a  center  rib  illustrating  the  failure  of  both  rib 
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a)  TYPICAL  RIB  FLANGE  FAILURE: 
SPECIMEN  BX-5B  (SKIN  REMOVED) 


b)  CENTER  RIB  ILLUSTRATING  FLANGE  FAILURE: 

SPECIMEN  BX-56 

FIGURE  75  TYPICAL  RIB  FLANGE  FATIGUE  FAILURES:  SPECIMEN  BX-5B 
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c)  RIB/FRAME  INTERSECTION:  SPECIMEN  BX-5B 

FIGURE  75  TYPICAL  RIB  FLANGE  FATIGUE  FAILURES: 
SPECIMEN  BX-5B  (CONCLUDED) 
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flanges.  This  failure  is  typical  of  designs  BX-4,  -5,  and  -6  in  that  rib  flange  fatigue 
was  experienced  on  flanges  next  to  the  front  and  rear  skins.  Figure  75c  is  a  closeup  of 
the  rib/frame  intersection.  No  fatigue  failures  of  the  rib  were  observed  at  locations 
other  than  the  flange  and  no  frame  fatigue  failures  were  observed  for  all  box  specimen 
structure  tested. 

The  obvious  conclusion  to  be  made  from  the  box  structure  fatigue  tests  is  the  importance 
of  designing  the  structure  so  that  the  rib  and  skin  are  highly  coupled  dynamically.  This 
will  be  achieved  by  selecting  rib  and  frame  thickness  identical  to  the  skin  thickness. 

To  emphasize  this  point.  Table  XIV  represents  a  summary  of  overall  rms  strain  response 
for  the  box  structure  designs  tested.  It  is  seen  that  the  strain  response  of  designs  BX-1, 
-2,  and  -3  is  more  uniformly  distributed  than  for  designs  BX-4,  -5,  and  -6. 

TABLE  XIV 

OVERALL  STRAIN  DISTRIBUTION  FOR  BOX  STRUCTURE 
(See  Figure  22) 


Overall  RMS  Strain,  Micro-inch/inch 
Strain  Gage  Number 


Specimen 

2 

A 

7 

10 

8 

JM 

9 

12 

J5 

BX-1  A 

580 

550 

— 

480 

500 

500 

500 

— 

— 

— 

BX-1B 

900 

610 

530 

— 

— 

— 

515 

— 

— 

— 

BX-2A 

495 

— 

530 

— 

— 

— 

— 

— 

— 

— 

BX-2B 

595 

560 

600 

— 

— 

490 

500 

— 

— 

— 

BX-3A 

420 

505 

500 

500 

470 

— 

430 

475 

— 

— 

BX-3B 

480 

495 

490 

— 

— 

140 

— 

— 

— 

490 

BX-4A 

520 

360 

580 

640 

680 

240 

350 

220 

320 

110 

BX-4B 

840 

600 

750 

1350 

1050 

450 

400 

720 

800 

500 

BX-5A 

— 

700 

1200 

2150 

2050 

— 

— 

1800 

1550 

800 

BX-5B 

— 

620 

800 

1200 

1050 

920 

— 

1500 

1300 

460 

BX-6A 

1150 

520 

520 

1500 

1100 

— 

700 

720 

580 

460 

BX-6B 

1050 

560 

950 

1800 

1750 

1600 

2000 

580 

IV  -  TEST  DATA  CORRELATION 


A.  Introduction 


The  analysis  presented  in  Section  II  suggests  ways  of  correlating  the  test  data  presented 
in  Section  III  for  the  stiffened  panel  stringer  specimens  and  the  nine-cell  box  structure 
specimens.  The  first  mode  response  was  in  general  the  most  predominant  mode  for  each 
type  of  structure.  For  the  stiffened  panel  specimens,  stringer  fatigue  failures  occurred 
at  the  clip  attachment  and  skin  fatigue  failures  occurred  at  the  frame  rivet  line.  Fa¬ 
tigue  failures  for  box  structure  occurred  at  the  frame  and  rib  rivet  lines  for  the  skin  and 
along  the  flange  rivet  line  for  the  ribs.  All  skin  fatigue  failures  for  both  types  of  struc¬ 
ture  fell  within  the  95%  confidence  limits  established  by  Ballentine  for  aluminum  plate 
fatigue  failure  (Reference  8,  Figure  34,  p.  91)  indicating  the  validity  of  Ballentine 's 
design  nomograph  (Reference  8,  Figure  73,  p.  151)  for  establishing  skin  thickness  for 
the  types  of  structure  considered  here . 

The  data  correlation  presented  here  is  concerned  with  stringer  fatigue  failure  at  the  clip 
attachment  and  flange  fatigue  failure  for  the  deep  plate-like  ribs  of  the  box  structure. 
The  approach  taken  for  data  correlation  is  somewnat  different  than  that  described  by 
Ballentine  in  that  overall  rms  stress  is  calculated  using  Equation  3  where  the  frequency 
and  static  stress  estimate  are  calculated  using  the  appropriate  expressions  from  the  anal¬ 
ysis  and  the  spectrum  level  of  the  excitation  and  damping  values  are  taken  from  the  ex¬ 
perimental  data.  The  resulting  values  of  overall  rms  stress  are  then  correlated  with  the 
measured  overall  rms  stress  (strain)  to  obtain  an  empirical  expression  for  overall  rms 
stress.  This  approach  is  open  to  question  since  scatter  exists  for  the  fatigue  data  (which 
is  assumed  to  be  error  free  for  this  correlation). 

B.  Stiffened  Panel  Stringer  Fatigue  Correlation 

Average  overall  rms  stress  levels  for  the  stringers  were  calculated  using  the  experimental 
data  presented  in  Table  X.  These  values,  denoted  by  a  e,  are  assumed  to  be  error  free. 
Values  for  the  response  frequency  of  the  (1,1)  mode  (stringer-torsion)  were  calculated 
using  Equations  25  and  26  and  values  for  the  stringer  stress  due  to  a  uniform  static  pres¬ 
sure  of  unit  magnitude  were  calculated  using  Equation  35  for  each  panel  design.  Using 
the  experimental  values  for  damping  (Table  IX)  and  rms  sound  pressure  spectrum  level 
(Table  X),  overall  rms  stress  levels  were  calculated  using  Equation  3.  These  values, 
denoted  by  a  c,  are  assumed  to  be  inaccurate  due  to  discrepancies  in  estimating  the 
frequency  and  static  stress  response.  Table  XV  lists  the  two  variables  a  c  and  ae  for 
the  various  specimen  designs.  These  data  are  plotted  in  Figure  76. 
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TABLE  XV 


> 


DATA  FOR  STRESS  CORRELATION:  STIFFENED  PANEL  STRINGERS 

Specimen 

g  x  10'2 
c 

G  X  10”2 
e 

STR-31 

4.582 

52.0 

STR-32 

8.420 

31.0 

STR-33 

1.942 

22.5 

STR-34 

4.597 

34.0 

STR-35 

1.498 

11.0 

STR-36 

3.820 

32.0 

STR-37 

23.330 

25.5 

STR-38 

6.050 

42.0 

STR-39 

22.450 

42.0 

STR-40 

28.440 

58.0 

STR-41 

6.790 

25.0 

STR-42 

2.900 

32.0 

From  the  data  listed  in  Table  XV,  the  correlation  coefficient  was  established  as  0.5754 
for  the  twelve  degrees-of- freedom.  This  level  of  data  correlation  is  not  very  significant 
(one  chance  out  of  fifty  that  there  is  no  correlation).  The  equation  for  the  regression 
line  is  established  as 

*=0.9(o)l/5,  ksi  (65) 

c  c 

The  regression  line  and  95%  confidence  limits  are  presented  in  Figure  76. 

C.  Box  Structure  Flange  Fatigue  Correlation 

Average  overall  rms  stress  levels  for  the  rib  flange  were  calculated  using  the  experi¬ 
mental  data  presented  in  Table  XIV.  These  values,  denoted  by  ~ae,  are  assumed  to  be 
error  free.  Values  for  the  frequency  of  Hie  (1,1,1)  mode  was  calculated  using  Equations 
44  and  45  and  the  static  stress  response  to  a  uniform  pressure  of  unit  magnitude  at  the 
mid  span  of  the  flange  was  calculated  using  Equations  48b  and  49.  With  these  data, 
modal  damping  data  (Table  XII)  and  rms  sound  pressure  spectrum  level  (T  ble  XIII),  the 
overall  rms  stress  levels  were  calculated  using  Equation  3.  These  values,  denoted  by 
7C,  are  assumed  to  be  inaccurate  due  to  discrepancies  in  estimating  the  frequency  and 
static  stress  response.  Table  XVI  lists  the  two  variables  ?c  and  <Te  for  the  various  spe¬ 
cimen  designs.  These  data  are  plotted  in  Figure  77. 
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FIGURE  77  CALCULATED  VERSUS  MEASURED  OVERALL  RMS 
STRESS:  BOX  STRUCTURE  RIB  FLANGES 
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TABLE  XVI 


DATA  FOR  STRESS  CORRELATION:  BOX  STRUCTURE  FLANGES 
(FLANGES  ON  THE  EXCITATION  SIDE) 


Specimen 

acx,° 

aex  10 

BX-IA 

2.2897 

5.974 

BX-1B 

2.2897 

7.365 

BX-2A 

1.5692 

5.279 

BX-2B 

1.5692 

6.154 

BX-3A 

2.8650 

4.738 

BX-3B 

2.8650 

5.996 

BX-4A 

13.5056 

6.625 

BX-4B 

15.2012 

11.580 

BX-5A 

20.2210 

21.130 

BX-5B 

22.7599 

20.300 

BX-6A 

24.8750 

13.583 

BX-6B 

24.8750 

16.840 

From  the  data  listed  in  Table  XVI,  the  correlation  coefficient  was  established  as  0.862 
for  the  twelve  degrees-of-freedom.  This  level  of  correlation  is  significant  (better  than 
one  chance  out  of  one  thousand  that  there  is  no  data  correlation).  The  equation  for  the 
regression  line  is  established  as 

ae=  0.245(ac)-411,  ksi  (66) 

The  regression  line  and  95%  confidence  limits  are  presented  in  Figure  77. 
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V  -  DEVELOPMENT  OF  DESIGN  EQUATIONS  AND  NOMOGRAPHS 


A.  Introduction 


The  information  required  to  produce  a  design  nomograph  for  a  specific  type  of  structure 
is  an  analytical  expression  for  the  stress  which  is  caused  by  the  random  acoustical  exci¬ 
tation  ana  the  allowable  fatigue  stress  and  life  relationship.  The  analytical  expressions 
developed  in  Section  II  for  frequency  and  stress  at  a  point  on  the  structure  were  modi¬ 
fied  empirically  by  the  data  correlation  in  Section  IV. 

To  be  useful  and  practical,  a  design  nomograph  also  must  be  as  free  as  possible  of  com¬ 
plexities,  easy  to  use,  and  give  an  accurate  representation  of  the  test  results  from  which 
it  was  derivea.  The  expressions  for  estimating  the  natural  frequency  of  the  fundamental 
mode  for  the  stiffened  flat  panel  and  box  structure  configuration  are  too  complex  to  be 
presented  in  the  form  of  a  nomograph.  Estimates  for  the  modal  stiffness  and  the  modal 
mass  can  be  readily  obtained  using  the  appropriate  equations  in  Section  II  for  any  mode 
desired.  For  extensive  parameter  studies,  computer  programs  have  been  developed, 
based  upon  the  theory  presented  in  Section  II,  for  estimating  the  frequencies  for  stif¬ 
fened  flat  panels,  box  structure,  and  a  three-cell  wedge  structure.  These  computer 
programs  and  a  description  of  the  data  input  format  are  presented  in  Appendix  IV. 

«  B.  Flat  Stiffened  Panels 


1 .  Skin  Design 

The  geometry  and  nomenclature  for  a  bay  of  stiffened  panel  structure  is  presented  in 
Figure  78.  All  skin  failures  were  observed  to  fall  within  the  confidence  limits  estab¬ 
lished  by  Ballentine  for  aluminum  skin  (Reference  8,  Figure  73,  p.  151).  Ballentine's 
design  nomograph  is  reproduced  as  Figure  79.  An  example  problem  is  included. 

2f  String?!-  Ftona?  §trs«« 

The  random  acoustic  pressures  on  the  surface  of  a  bay  of  stiffened  panel  structure  are 
transferred  to  the  supporting  structure  predominantly  by  a  transverse  shear  loading.  This 
loading  causes  the  tnin-walled  open-section  stiffeners  to  bend  and  twist.  The  stringer 
loading  is  reacted  along  the  skin-stringer  rivet  line  and  the  clip  attachments  of  the  frame 
and  stringer.  The  stringer  stress  AT  THE  CLIP  ATTACHMENT  is  estimated  by  the  equation 


-J{U/  C  ,)hb3 

7c  ”  i*  [(bA)  +  (aA))  Vfi> 


(67) 


*  =  (I 

XX  zz 


zz 


for  the  fundamental  mode,  f],  of  the  panel.  The  frequency  of  the  fundamental  mode, 
fj,  is  obtained  by  using  Equations  22  and  23  for  m  =  n  =  1.  The  damping  ratio,  C  ] , 
can  be  estimated  as  0.01  to  0.02  in  the  absence  of  experimental  data. 
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The  sound  pressure  spectrum  level,  S  (fj),  for  the  fundamental  mode  is  expressed  in  units 
of  pressure  as  K 

S^)  =  2.91  x  10<SPl/20-9>  ,  lbs/in2/ -/  Hz  (68) 

where  spl  is  the  spectrum  level  in  dB 

From  the  empirically  derived  expression  for  stress.  Equation  65,  the  estimated  stress  in 
the  stringer  flange  at  the  clip  attachment  is  obtained. 

EXAMPLE:  Consider  a  stiffened  panel  with  the  following  data  - 

a  =  9.0  in  b  =  18.0  in  t  =  0.040 

h  =  1.25  in  I*  =  0.01255  inf  (-4  zee  section.  Table  VI) 

fj  =  148  Hz  (from  Equations  22  and  23) 

Sp  =  0.01158  (from  Equation  68,  spl  =  132  dB.) 

Cj  =  0.015  (estimated) 

(b/a)+(a/b)  =  2.50 

Then,  ~c  =  5.16  (from  Equation  67) 

and  ae  =  4.64  ksi  (from  Equation  65). 

From,  Figure  44,  the  cycles  to  failure,  N-,  estimated  to  be  1.22  x  10^.  From  Figure 
79,  and  the  above  data  the  number  of  cycfes  to  failure  are  estimated  to  be  2.2  x  10' 
cycles. 

The  above  procedure  should  be  repeated  until  the  desired  fatigue  life  for  the  skin  and 
the  stringers  is  achieved. 


To  calculate  the  stresses  in  the  skin  and  rib  flanges  for  a  box  structure,  the  design  pro¬ 
cedure  utilizes  Timoshenko's  approximate  technique  for  computing  stresses  in  continuous 
supported  panels  (Section  II. D. 4),  the  frequency  estimate  for  the  (1,1,1)  mode  (Equa¬ 
tions  44  and  45),  and  the  empirically  determined  expression  for  stress  given  by  Equation 
66.  The  method  is  best  illustrated  by  an  example. 

EXAMPLE:  Consider  the  data  for  specimen  BX-1  (Table  VII).  From  Equations  48a  and 
48b,  the  bending  moment  at  the  middle  of  the  short  side  of  the  center  bay  is  given  by 

My  =  ((0. 0575)(1 00)  +  (0.  Q575)(l  00))/3  =  3.84 


/ 


and  at  the  middle  of  the  long  side  of  the  center  bay 

Mx  =  ((0.065)(121)+ (0.082)(1 00))/3  =  5.35 

for  a  uniform  pressure  of  unit  magnitude  over  the  surface  of  the  structure. 

To  obtain  stress  in  the  skin,  frame  flange,  and  rib  flange,  one  uses  Equation  49.  For 
the  skin  at  the  middle  of  the  short  side  (frame  line)  of  the  center  bay  the  skin  stress 
response  is 

a,  =  (6.0)(3.84)/(0.063)2  =  5800  psi, 

and  the  frame  flange  stress  response  is 

°2  =  (6 . 0)(3 . 84)/(0 . 080)2  =  3600  psi. 

For  the  skin  at  the  middle  of  the  long  side  (rib  line)  of  the  center  bay  the  skin  stress 
response  is 

°3  -  (6 . 0)(5 . 35)/(0 . 063 )2  =  8080  psi, 
and  the  rib  flange  stress  response  is 

=  (6 . 0)(5 . 35)/ (0 . 063 )2  =  8080  psi. 

From  Equations  44  and  45  the  frequency  of  the  fundamental ,  (1,1,1),  mode  is  estimated 
to  be  79.0  Hz.  For  this  example,  the  damping  ratio  is  estimated  to  be  0.020  for  the 
entire  structure,  and  a  sound  pressure  spectrum  level  of  124  dB/VHz.is  assumed.  From 
Equation  3,  one  obtains  the  expression  (ignoring  the  constant) 

*ci  ■  «w/ci>  yfi>v  <69> 

Then,  for  the  above  data 

f,  =  79.0,  Hz.  C,  =  0.020  Sp(f t )  =  0.00462 

so  that 

*  j  =  (62. 8)(. 00462 )(5800)  =  1682  psi, 
ac2  =  (62. 8)(. 00462)(3600)  =  1045  psi, 
oc3  =  (62. 8)(. 00462 )(8080)  =  2340  psi, 
ac4  =  (62. 8)(. 00462X8080)  =  2340  psi. 

From  Equotipq  66,  the  estimated  overall  rms  stress  is  given  by 
ael  =  0.245  (1682)*41 1  =  5.20  ksi, 

■ae2  =  0.245  (1045)'411  =  4.25  ksi, 

ae3  =  0.245  (2340)*411  =  6.30  ksi, 

a#4  =  0.245  (2340) ,411  =  6.30  ksi. 


122 


For  the  skin  stresses  o  ^  and  aeq  the  number  of  cycles  to  failure  are  obtained  from 
Figure  80  as  1.2  x  10'  cycles  ana  5  x  10^  cycles,  respectively.  For  the  frame  and  rib 
stresses  the  number  of  cycles  to  failure  are  obtained  from  Figure  69  as  2.3  x  10?  cycles, 
respectively. 

The  above  method  can  be  repeated  until  the  desired  fatigue  life  is  obtained. 


D.  Frequencies  of  a  Three-Cell  Wedae  Structure 


The  vibration  analysis  of  a  three-cell  wedge  structure  is  presented  in  Section  II. E.  An 
example  will  be  presented  here  for  a  specific  data  case.  The  results  given  here  were 
obtained  using  the  computer  program  given  in  Appendix  IV  (WDGVIB);  however.  Equa¬ 
tions  57  and  58  can  be  used  if  desired. 

EXAMPLE:  With  reference  to  Figure  9,  consider  the  data  case: 
jts  =  2*s  =  t  =  0.050  in. 
a  =  bj  =  ^2  =  ^3  =  in. 

0-  =  11 .46  degrees  (0.2  radian) 

Y  =  0.101  lbs/in2 
E  =  1.03  x  107  psi 
v  =  0.302 

]Dj  =  114.46  in  -  lb. 

Then,  the  first  few  modal  frequencies  are 


mode  number 
(m. n,q) 

f  ,  H 
mnq' 

(M/0 

134.0 

(M/3) 

298.0 

0/2,1) 

298.0 

(1,2,3) 

611.0 

0,3,1) 

530.0 

0,3,3) 

964.0 

(2,1,1) 

261.0 

Typical  rib  mode  shapes  are  illustrated  in  Figure  11  for  this  data. 
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FIGURE  80  SKIN-STRINGER  PLATING  FATIGUE  CURVE 
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E.  Curved  Panel  Frequencies 


1.  Szechenyi's  Approximate  Formulae 

Application  of  Szechenyi's  approximate  formulae  for  estimating  the  natural  frequencies 
of  cylindrical  panels  will  be  illustrated  by  an  example.  The  data  used  is  obtained  from 
experimental  values  quoted  by  Petyt,  Reference  13. 

EXAMPLE:  Using  the  frequency  expression  giver,  by  Equation  59,  and  Figures  13  and  14 
one  can  compute  frequencies  for  the  various  lower  order  panel  modes  of  a  cylindrical 
panel.  Consider  the  data  case  (see  Figure  12) 


a  =  3.0  in,  b  =  4.0  in,  t  =  0.013  in,  R  =  30.0  in 

E  =  107  psi,  Y  =  0.096  lb/in3,  v«  0.33 

The  comparison  between  Szechenyi's  approximate  formula  (Equation  59)  and  the  experi¬ 
mental  results  quoted  by  Petyt  are 

Equation  59  Experiment 


fll  ■ 

1053., 

Hz. 

fl2  “ 

968., 

Hz. 

f13  = 

1192., 

• 

N 

X 

f.4  ■ 

1730., 

Hz. 

—  ,  Hz. 

814.,  Hz. 

940.,  Hz. 
1735/1770,  Hz. 


Considering  the  simplicity  of  Szechenyi's  method  the  results  are  considered  to  be  good. 

Zt-PJynTblgg's  Approximate  Pgrnmiqg 

The  expressions  for  the  frequency  ratio  and  stress  ratio  for  curved  to  flat  panel  response 
in  the  fundamental  mode  presented  in  Section  II. F. 2  can  be  applied  directly  for  specific 
data  cases.  The  nomogram  developed  by  Plumblee  (Figure  87,  Reference  8,  page  174) 
is  presented  as  Figure  81 .  This  nomogram  is  used  in  conjunction  with  Figure  79  for  esti¬ 
mating  curved  panel  fatigue  life  for  acoustic  oxcitation. 


FIGURE  81  EFFECT  OF  CURVATURE  ON  STRESS  RATIO  FOR  SKIN-STRINGER  PANELS  WITH  RANDOM  LOADING 
(  REFERENCE  8  ) 


IV  -  CONCLUSIONS 


Methods  have  been  developed  to  estimate  the  acoustic  fatigue  resistance 
of  aircraft  structural  components.  These  methods  apply  to  flat  stiffened 
panels,  and  to  box  type  and  wedge-type  structures  such  as  flaps  and  trailing 
edges  of  wings  and  control  surfaces.  In  the  acoustic  fatigue  analysis  the 
response  frequencies  of  these  structural  components  are  estimated  by  the  use 
of  the  Rayleigh  energy  method  which  accounts  for  the  bending  and  torsional 
flexibilities  of  the  supporting  structure.  The  mean  square  stress  response 
of  the  skin  which  is  subjected  to  the  random  acoustic  load  is  computed  by 
the  use  of  Miles  equation  for  a  s i ngle-degree-of-f reedom  system.  The  life 
of  the  structure  is  then  estimated  by  using  experimentally  derived  S-N 
curves  and  a  Miner's  cumulative  damage  calculation.  Acoustic  fatigue  tests 
of  twelve  stiffened  flat  panels  and  six  box  structure  designs  were  performed 
to  check  the  theory  and  supply  the  S-N  data  for  the  analysis.  The  resulting 
S-N  curves  are  given  in  two  graphs,  one  for  the  stiffened  plates,  the  other 
for  the  box  structure.  Computer  programs  are  also  given  for  estimating 
frequencies  of  the  stiffened  flat  panels,  box  structures,  and  wedge 
structures.  Based  upon  this  program  the  following  conclusions  were  made: 

a.  The  results  of  the  analytical  and  experimental  investigations 
reported  here  indicate  that  the  overall  rms  stress  induced  in  stiffened 
panel  skins  and  box  structure  skin  and  skin-rib  connections  by  broad-band 
noise  can  be  estimated  with  sufficient  accuracy  for  design  purposes  using 
approximate  results  from  thin  plate  theory  and  Miles'  response  theory. 

b.  The  approach  taken  to  establish  stringer  flange  design  data  is 
considered  to  be  valid  within  the  nature  of  the  data  obtained  in  the  experi¬ 
mental  program.  The  difficulties  in  measuring  strain  at  the  failure  point 
and  the  extreme  strain  gradients  present  at  the  clip  attachment  produce  an 
additional  scatter  to  the  fatigue  data  which  is  inherent  in  the  type  of 
structure  considered. 

c.  The  confidence  limits  established  by  the  limited  fatigue  data 
obtained  in  the  experimental  program  are  too  broad.  Additional  testing 
will  be  required  to  better  define  the  acoustic  fatigue  characteristics  for 
both  stringer  flange  failures  and  rib  flange  failures. 

d.  The  vibration  analysis  for  wedge  structure  should  be  reformulated 
to  provide  both  frequency  and  stress  response  estimates.  Even  through  the 
frequency  analysis  presented  here  yields  reasonable  results,  the  simple 
mode  shapes  used  do  not  allow  an  estimate  of  stress  to  be  made.  Acoustic 
fatigue  testing  of  wedge-type  structure  will  be  required  to  establish 
fatigue  design  criteria. 

e.  The  design  equations  developed  from  the  analytical  and  experimental 
programs  will  increase  acoustic  fatigue  design  capabilities,  but  these 
equations  should  be  considered  preliminary  until  more  experimental  data  is 
avai lable. 
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APPENDIX  I 

GENERAL  EXPRESSIONS  FOR  STIFFENER  GEOMETRIC  PARAMETERS 


The  geometric  parameters  defined  here  are  developed  in  terms  of  a  centroidal  (x,y,z) 
coordinate  system.  General  expressions  for  the  cross-sectional  area,  the  area  moments, 
torsion  constant,  and  warping  constants  are  presented  for  zee,  channel,  and  hat 
cross-section  shapes.  The  definition  of  these  parameters  is  as  follows: 


x,  the  location  of  the  centroid,  as  indicated 
e,  the  location  of  shear  center,  as  indicated 
A,  the  cross-sectional  area 


1  =  f  z2dA 

XX  J 

1  -  f  xzdA 

xz  J 

A 

A 

j  -  /[(<p,x  -  *)2 + 

(cp,z  +  x)2  ]  dA 

A 

REz  =  /  XCpdA 

REx  =  /  ZCpdA 

A 

A 

A 

(St.  Venant's  Torsion  Constant) 

v  h2“ 

A 


where  cp(x,z)  is  the  warping  function  for  the  cross-section  with  the  pole  taken  as  the 
shear  center . 


* 

4 


i 
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GENERAL  EXPRESSIONS  FOR  STIFFENER  GEOMETRIC  PARAMETERS 


A  =  t(h  +  2b) 

,  s  J_  [h2(6b  +  h)  +  t2(3h  +  2b)] 

xx  12 

|  =-t(2b+t)(2b-t) 

XZ  o 

lzz  =  ^[8b3+ht2] 

t3 

j=L[2b+h] 

„  tb3h2(b  +  2h) 
e  ”  12(2b+hT 


FIGURE  A-l-1  GEOMETRIC  PROPERTIES  -  ZEE  SECTION 
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GENERAL  EXPRESSIONS  FOR  STIFFENER  GEOMETRIC  PARAMETERS 


x  =  b2/(2b  +  h) 

e  =  3b2/(6b  +  h) 

A  =  t(h  +  2b) 

h^6b  +  +  +  2b^ 

I  =0 
xz 

I  =  [  12h^2  +  8b3  -  24xb(b  -  x) 

J  =  i  t3(2b  +  h) 

_  »b3h2(3b  +  2h) 
e  1 2(6b  +  h) 


hbH 

H- 

i-  1 

1 

+ 

E  t* 

1  X  1 

-A- - ►  h 

+■  1 

L 

12b(b  -  2x)t  +  6(b  -  x)t2  +  t3] 

NOTE:  For  warping  constants,  the  pole  is  taken  at  the  shear  center. 


FIGURE  A-l-2  GEOMETRIC  PROPERTIES  -  CHANNEL  SECTION 
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GENERAL  EXPRESSIONS  FOR  STIFFENER  GEOMETRIC  PARAMETERS 


b(b  +  2hj) 

*  =  +  K  +  2f^ 


b(3bh2  +  2h,(3h2  -  4hft) 

[h2(6b  +  h)  +  2h  ^  (3h^  +  6hh1  +  4h^)] 


e  - 


A  =  t(h  +  2b  +  2h,) 

.  '  t(h  -  „3  4  Kb  ♦  DO2  +  3b2)  -  j|  «*,  *  '>  «»l  * ,)2  +  ’2(h  +  ^ 

xx  12  0 


=  0 


xz 


4  kk  -  .H.2  *  §  Kb + «<b  ♦  >k- w 


zz 


+  lt(2h1  -t)[t2+24(b-^)2] 


j=^t3[h+2b+2h1] 


R  =0 
ex 


R  =0 
ez 


r  =J^  [2bV  +  6bh2e(e  -  b)  +  2b2h1(3h2  -  6HH,  +  4h] 


-  4bh1«(3h2  -  4h^)  +  e2(h  +  2h1)3] 


FIGURE  A-l-3  GEOMETRIC  PROPERTIES  -  HAT  SECTION 
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APPENDIX  II 


TEST  DETAILS  AND  INSTRUMENTATION 


A.  Introduction 

In  general;  the  test  details  and  instrumentation  used  to  conduct  the  experimental  program 
are,  except  for  the  statistical  analysis,  identical  to  the  procedures  described  by  Ballentine 
(Appendix  I,  Reference  8).  The  procedures  used  for  the  box  specimen  mode  studies  were 
slightly  different  than  that  used  for  the  stiffened  panel  specimens  due  to  the  box  structure 
configuration.  These  differences  are  described  in  Section  III.  The  statistical  data  analysis 
will  be  described  here,  and  the  reader  can  obtain  other  details  as  described  in  Reference  8. 


Statistical  Data  Analysis 


The  data  used  to  obtain  the  S-N  curves  presented  in  this  report  consisted  of  R.M.S.  values 
for  strain,  zero  crossings  with  negative  slope  of  the  strain  signal,  and  the  time  that 
fatigue  failure  was  first  observed.  All  fatigue  failure  data  reported  in  this  report  corre¬ 
spond  to  a  failure  initiated  at  oc  very  close  to  the  strain  gage.  The  zero  crossing  count 
was  obtained  using  a  B  &  K  Model  161  Probability  Density  Analyzer,  a  frequency  counter, 
and  the  recorded  strain  signal.  The  average  of  the  frequency  count  (zero  crossings  with 
negative  slope)  for  three  minutes  of  strain  data  was  taken  as  the  zero  crossing  count.  The 
block  diagram  for  the  system  is  illustrated  in  Figure  A-ll-1 . 


The  amplitude  distribution  plots  illustrated  in  this  report  were  also  obtained  using  the 
B  &  K  Model  161  Probability  Density  Analyzer.  Six  minutes  of  recorded  strain  data  were 
used  to  obtain  each  plot.  The  Gaussian  distribution  referenced  in  each  plot  was  obtained 
by  analyzing  a  1  volt  rms  random  noise  signal  from  a  B  &  K  Model  1024  sine-random  noise 
generator.  The  block  diagram  for  the  system  is  illustrated  in  Figure  A-ll-2. 
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APPENDIX  III 


EXPERIMENTAL  SPECIMEN  MODE  SHAPES 


I .  Introduction 

The  experimental  frequencies  and  mode  shapes  described  here  were  obtained  as  described 
in  Sections  III.C  and  III.D.  The  specimen  plan  view  is  illustrated  (no  scale).  The  frame 
and  stringer  attach  lines  are  indicated.  The  frame  lines  are  horizontal  and  the  stringer 
attach  lines  are  vertical.  The  node  lines  for  a  given  bay  are  indicated  by  dashed  lines. 
Predominant  modes  for  each  specimen  are  indicated  by  an  asterisk. 
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FIGURE  AIM-  STIFFENED  PANEL  MODE  SHAPES 
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FIGURE  AMI-  STIFFENED  PANEL  MODE  SHAPES 
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FIGURE  AIM-  STIFFENED  PANEL  MODE  SHAPES 
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STIFFENED  PANEL  MODE 
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FIGURE  Alll-  STIFFENED  PANEL  MODE  SHAPES 
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FIGURE  AIII-  STIFFENED  PANEL  MODE  SHAPES 
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FIGURE  AIII-  STIFFENED  PANEL  MODE  SHAPES 
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APPENDIX  IV 


COMPUTER  PROGRAMS 


1 .  Introduction 


The  computer  programs  listed  here  are  based  upon  the  analysis  given  in  Section  II.  All 
programs  compute  frequencies  and  modal  amplitudes.  The  wedge  structure  vibration 
analysis  computes  the  mode  shapes  for  the  ribs  as  a  function  of  position  along  the  radial 
axis  at  the  center  line  of  the  rib.  These  computer  programs  have  been  written  for  the 
UNIVAC  418  computer;  however,  with  very  slight  modifications,  they  can  be  adapted 
to  any  computer.  Special  functions  and  subprograms  are  listed  with  the  main  program. 

2.  Nine  Bay  Stiffened  Panel  Vibration  Analysis  (PLTVIB) 

Purpose:  To  compute  the  natural  frequencies  and  modal  amplitudes  for  the  panel  bays 
and  ribs  as  described  in  Section  II. C.  The  program  is  listed  in  Figure  A  I V- 1 .  Typical 
output  format  is  listed  in  Figure  A  IV —2 . 

Subprograms  Required:  SGN  (K) 

Input  Data  Format  (PLTVIB):  5  Cards  per  Data  Case 


CARD  1 


COL  (FORMAT) 

1(14) 

5(12) 

7(12) 

9(12) 

5(12) 

NAME 

NDATA 

MS 

MF 

NS 

NF 

CARD  2 


COL  (FORMAT) 

1 (El  0.4) 

11 (El  0.4) 

21 (E10.4) 

31 (El  0.4) 

NAME 

A1 

A2 

B1 

B2 

CARD  3 


COL  (FORMAT) 

1  (El  0.4) 

1 1 (El  0-4) 

EQ2QQ1 

NAME 

TS 

GMS 

ES 

PRS 

CARD  4/5  (Cont'd) 


COL  (FORMAT) 

1  (E10.4) 

1 1  (El  0.4) 

21  (El  0.4) 

31  (El  0.4) 

NAME 

gx/gy 

XJ/YJ 

EX/EY 

GMX/GMY 
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CARD  4/5 


COL  (FORMAT) 

41  (E10.4) 

51 (El  0.4) 

61  (E10.4) 

NAME 

AX/AY 

XI/YI 

WCX/WCY  ! 

NDATA 

MS,MF 

NS,NF 

A1,A2 

B1,B2 

TS 

GMS 

ES 

PRS 

GX/GY 

XJ/YJ 

EX/EY 

GMX/GMY 

AX/AY 

XI/YI 

WCX/WCY 


four  digit  data  identification  number 
initial  and  final  values  for  mode  number  m 
initial  and  final  values  for  mode  number  n 
panel  bay  dimensions  in  X-direction  (see  Figure  1) 
panel  bay  dimensions  in  Y-direction  (see  Figure  1) 
skin  thickness 

weight  density  of  skin  material 

Young's  modulus  of  skin  material 

Poisson's  ratio  of  skin  material 

shear  modulus  of  stringers  parallel  to  X/Y  axis 

St.  Venant's  torsion  constant  for  stringers  parallel  to  X/Y  axis 

Young's  modulus  for  stringers  parallel  to  X/Y  axis 

weight  density  of  material  for  stringers  parallel  to  X/Y  axis 

cross  section  area  for  stringers  parallel  to  X/Y  axis 

Ip^/A^  and  l^^/A  (see  Equations  8a  and  8b) 

warping  constant  ^or  cross-section  of  stringers  parallel  to  X/Y  axis. 


3.  Nine-Cell  Box  Structure  Vibration  Analysis  (BOXVIB) 

Purpose:  To  compute  the  natural  frequencies  and  modal  amplitudes  for  the  skin,  rib,  and 
frame  of  a  nine-cell  box  structure  as  described  in  Section  II. D.  The  program  is  listed  in 
Figure  A  IV-3  and  typical  output  is  listed  in  Figure  A  IV-4. 

Subprograms  Required:  STIFF(SK,  I,  J,K),  MASS(SM,I,  J,K),  F(N,I,J),  SGN(K). 

Input  Data  Format  (BOXVIB):  Five  Cards  per  Data  Case 


CARD  1 


COL  (FORMAT) 

1(14) 

5(12)  I 

NAME 

NDATA 

IOUT 

CARD  2 


COL  (FORMAT) 

KI2) 

3(12) 

5(12) 

6(12) 

7(12) 

11(12) 

NAME 

MS 

MF 

NS 

NF 

QS 

QF 
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CARD  3 


COL  (FORMAT) 

1  (El  1.4) 

12(E1 1 .4) 

23(E)  1.4) 

34(E1 1 .4) 

45(E 1 1 .4) 

NAME 

Al 

A2 

B1 

B2 

H 

CARD  4 


COL  (FORMAT) 

1  (El  1.4) 

12(E1 1  .A) 

23(E1 1 .4) 

34(E1 1 .4) 

45(E1 1 .4) 

NAME 

T1S 

T2S 

GMS 

ES 

PRS 

CARD  5 


COL  (FORMAT) 

1 (El  1 .4) 

12(E1 1 .4) 

23(E1 1 .4) 

34(E1 1 .4) 

45(E1 1 .4) 

NAME 

T1R 

T2R 

GMR 

ER 

PRR 

NDATA 

IOUT 

MS,MF 
NS,NF 
QS,QF 
Al ,  A2 
B1 ,82 
H 

T1S/T2S 

GMS 

ES 

PRS 

T1R/T2R 

GMR 

ER 

PRR 


four  digit  data  identification  number 

output  option:  =  0,  frequencies  are  printed  for  each  mode 

=  1,  frequencies  and  modal  amplitudes  are  printed  for 
each  mode 

initial  and  final  values  for  mode  number  m 
initial  and  final  values  for  mo  3  number  n 
initial  and  final  values  for  mode  number  q 
length  of  panel  bays  parallel  to  x-axis  (see  Figure  3) 
length  of  panel  bays  parallel  to  y-axis  (see  Figure  3) 
height  of  box 

upper/lower  skin  thickness  (see  Figure  3) 
weight  density  of  skin  material 
Young's  modulus  for  skin  material 
Poisson's  ratio  for  skin  material 
rib/frame  thickness  (see  Figure  3) 
weight  density  of  rib  and  frame  material 
Young's  modulus  of  rib  and  frame  material 
Poisson's  ratio  for  rib  and  frame  material 


4.  Three-Cell  Wedge  Structure  Vibration  Analysis  (WDGVIB) 

Purpose:  To  compute  the  natural  frequencies  and  modal  amplitudes  for  the  ribs  and  skin 
of  a  three-cell  wedge  structure  as  described  in  Section  II. E.  The  computer  program  is 
listed  in  Figure  A  !V-5and  typical  output  is  listed  in  Figure  A  IV-6. 
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Subprograms  Required:  SGN(K) 

Input  Data  Format  (WDGVIB):  5  Cards  per  Data  Case 


CARD  1 


COL  (FORMAT) 

1(14) 

5(12) 

7(12) 

9(12) 

11(12) 

13(12) 

15(12) 

NAME 

NDATA 

MS 

MF 

NS 

NF 

QS 

QF 

CARD  2 


COL  (FORMAT) 

1  (El  0.4) 

1 1  (El  0.4) 

21  (El  0.4) 

31 (El  0.4) 

NAME 

A1 

B1 

B2 

ALPHA 

CARD  3/4 


COL  (FORMAT) 

1  (El  0.4) 

11 (El  0.4) 

2 1  (El  0.4) 

3 1  (El  0.4) 

NAME 

T1S/T2S 

GM1S/GM2S 

E1S/E2S 

PR1S/PR2S 

CARDS 


COL  (FORMAT) 

1(E10.4) 

1 1 (El  0.4) 

2 1  (El  0.4) 

31 (El  0.4) 

NAME 

T1R 

GMR 

ER 

PRR 

NDATA 

four  digit  data  identification  number 

MS,MF 

initial  and  final  values  for  mode  number  m 

NS,  NF 

initial  and  final  values  for  mode  number  n 

QS,QF 

initial  and  final  values  for  mode  number  q 

NOTE:  only  odd  values  of  q  are  considered.  QF  ^  5. 

A1 

radius  of  rib 

B1/B2 

rib  spacing  (B3  =  Bl:  see  Figure  9) 

ALPHA 

TOTAL  wedge  angle  in  DEGREES 

T1SA2S 

thickness  of  upper/lower  cover  sheet  (see  Figure  9) 

GM1S/GM2S 

weight  density  of  upper/lower  cover  sheet  material 

E1S/E2S 

Young's  modulus  of  upper/lower  cover  sheet  material 

PR1S/PR2S 

Poisson's  ratio  of  upper/lower  cover  sheet  material 

T1R 

thickness  of  rib 

GMR 

weight  density  of  rib  material 

ER 

Young's  modulus  of  rib  material 

PRR 

Poisson's  ratio  of  rib  material 

166 


Tulb  PROGRAM  EV  LuATEs  the.  modal  mass  and  stiffness, 
FhEQUENCY  and  mode  shapes  of  a  STIFFENED  panel 
Dimension  a<?3) >F(5> ,*k(3) #yk(3> #w(3#t) #Sx(2#3) #sy(2#3) 
In=1 
Iq=2 

REAO(lN»100)NDArA,MS»^F»NS»NF 

KeAD(IN.101)A1*A2,Bi,b2 

Rfc.ADllN,10l)TS»GMS»ES»PRS 

read  < in# ioi i  gxj.*  J  #  Ex#  ^x»_a.Xi_xjL#  wc.x _ 

KlAD(  lN»iol)GY»  YJ»Ey#GI*iY  #aY#YI#WCy 
A(1)=A1/A2 
A(2)=Bl/0£ 

A(3)=A2/02 
A(4)=A2*B2 

A(5)=A(4)**3/B2  _ _ _  _ _ 

A(6)=A(5)*B2/A2 
A(7)=A(1)**3 
A ( 8) -A ( 2 ) **3 

A(9)=( A2*2.*Al)*(b2+2.*Bl)*GM5*TS/4. 
A«10)=9.ttb959*GMX*Ay*xI/B2*A(3)*(l*+2.*A(7> ) 

A{il)  =  9.8b959*GHY*AmX/A2'A<3)*<l«^2..*A<a>-L_ 

A (12 ) =ES*  TS**3/ ( 1 ? . « ( j . -PRS*PRS )) /A ( 4 ) *2U . 35227 

A(13)=EX*«Cx/A(5!*961.3B5u6 
A(l4)=tY*wCY/A(t»)  *961.3faB(i6 
A(15)=61/A1 
A'lb)=ril/A2 
A(17>=u2/A1 
A(la)=GX*xJ/(EX*WGX) 

A(19)=GY*VJ/(EY*WCY) 

A(20)=A1**2 
A(2l)=A2**2 
A(22)=B1**2 
A(23)=c2**2 
DO  202  MrMSiMF 
SflMsSGMM) 

F-t=FLOAT(M) 

Ffi  2=FM*FM 
Fn>4=F,V2*FM2 
P1MA2=3.14I59*FM/a2 
F;  PlSQ=(FM*3,l4l5q)  **; 

Du  201  N=nS,NF 
St,N=SGN(N) 

Ff.=FuOAT(N) 

Fn2=Fn*FN 
Fr14=F(,2*FN2 

PlNa2=3.l4l59*FW/b2 

F  OFNsFM/FN 


FIGURE  A-IV-1  PROGRAM  PLTVIB 


^ ( i )  =  (  (FM0Fn*A(i5)M./(A<i5)*FM0Fn)  )*e2)*4»*A(l)*A*2) 

F(2)  =  (  (FM0FN/ACl)4ALaiyFHQFN>)**2  _ 

F<3)=< (FM0FN*A(l6)+l./(A(i6)*FMbFN) )**2)*2»*A(2) 

F (4 )  =  ((FMOFN»A(i7)»t./(A(l7)»FM0Fh!)  )»»2)»P«»A(1 ) _ 

F  <  5 )  =F  ( 1 )  *F  ( 2 )  «-F  ( 3 )  +F  { 4 ) 

FkPISQ=(Fn*3.1415c))**2 

Xk<1)=A(18)*A(2U) 

Xk<2)=A(18)*A(21) 

XK(3)  =  l.^XK{2)/FMpISQ+2./A(l)*(!..^XK(1  )/FmPISQ) 

Yk ( 1 ) “A { 19) *A (22 ) _ _ 

Yk(2)=A(19>*A(23) 

YK(3)=l.+YK(2)/FNpIS0^2./A(2)*(i*+YK(i)/FNPlSQ) 

Sf-  M=  ( A  ( 9)  +A  ( 10 )  *Fm2+ A  { 1 1 )  *Fm2  )  /386  • 

SkM=A { 12 ) *FM2*FN2*F ( 5 ) +A ( l 3 ) *FM4*FN2*xK ( 3 ) ♦A ( 14 ) *FM2*FN4*yK ( 3 ) 
FhEQ=SGRT ( SKM/SMM )  /ft .  ?fl 31fl 
.  wa#U=SNM*SNN*A(l)»Ai2) 

W  ( 2 » 1 ) =SNN* A ( 2 ) 

W ( 3i 1 ) =W ( 1 *i) 

W(1»2)=SNW*A(1) 

W (2? 2) =1  • 

*(3»2)=W(1#2) 

n  ( 1 1 3)  ( 1  r  l  J  „  _  _  .  ...  . 

W(2»3)=W<2»1) 

W(3i3)=W(1'1) 

Sx ( 1 » 1 ) =SNM*PlNti2*A { 1 ) 

Sx(1»2)=P1NQ2 

Sx<l*3)=Sx(l»l) 

Sx ( 2 » 1 ) =5nN*SX (  IjjJ . .  .  .  . . 

Sx(2»2)=SNN*PlNB2 

Sx(2»3)=Sx(2»l) 

Sy ( 1 # 1 ) =SNN*PIMA2*A { 2  j 

Sy(1»2)=PIMA2 

Sy(1»3)=SY(1»1) 

SYUfi)=SNM*SYU»iJL _  _ _ 

Sy(2»2)=SNM*PIMA2 

Sy(2»3)=SY(2»l) 

*rITE(IO»215) 

mrITE(I0>22U)N0ATa 

«(hITEU0»225)M»n»FREG 

WRITE  ( 10 » 2301  SMM»-S&4 _  _  _  . . . . 

WrITE( 10*235) 

write (io»24a) ( (w(i»j),j=i»3) » 1=1*3) 

WRlT£(lOf245) 

Write (I0»240) l 'SX{I,J) ,J=i»3>»I=1»2) 

W*<ITE(IO»250) 

WhITe(I0»240M(SYa«JliJ=U,l)iI=If2)  . .  .... 

?oi  Continue 


FIGURE  A-IV-1  PROGRAM  PLTVIB  (Continued) 
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GC  TO  200 
ino  F(,RMAT(I4,4I2> 

101  FrRMAT(7LlO,4) 

?lb  FuRMT(lHl»9X»42H,>Y!^ic  aNmLYSIS  OF  /  STIFFENED  FLAT  PANEL) 

220  Fl,H;;AT<//.23X.9HDaTa  CASE,l5) 

>25  FvjKF-AT  (//»3a»13HMqDC  uUMbER  <  #  12#  lH, ,  j2 » 1H> » 10X  #  1qHFREqUEnCY=.E.12, 
lb»4milZ* ) 

?3li  FcHsAT  (/»3XillHrt00AL  aSS=#L12 . 5# 3X » KHMOpAL  STIFfMESS=»Ei2,5) 

^35  F-RM1  (//#22X»lbH,iOr  al  AMpLi TUDES , ///? 3X » i 4hPANEL  -  W(I»J)»//) 

?4C  Fo^Ar(10x»£l2.b»,iX#(_i2.5,2x#El2.5> 

°4f>  FcRNAT(//»l0X»39M‘iTlFf  LnErS  IN  THE  X  DIRECTION  -  s*(I»J>#//> 

■>50  F  qR.'IAT  ( // » 10X » 39HSTIF  FLHLkS  IN  THE  Y  [.IREcTlON  -  sY(I#J)#//) 

Si  OP 

fc;.J 


Function  sgnin) 

Csl.O 

Do  5  1-1 »K 
Cr-C 

b  CcNTIN'JE 
S(,N=C 
Rj-  TURN 
Ci.D 


FIGURE  A-IV-1  PROGRAM  PLTVIB  (Concluded) 


dynamic  analysis  of  a  stiffened  flat  Panel 


data  case  37 

VODE  NUMBER  (  1»  1)  FREQUENCY:  0.23672E+03»H^. 

rO^AL  MASS:  Q.16264E-02  MODaL  STIFFnESS=  0.35979E+04 

MOjAl  amplitudes 
PaNfl  -  W(I»U) 


0.20370E+00  -C.faiUlL+UO  0.20370E+00 
-0.33333E+G0  0 .  ltfiOOE+Ol  -0.33333E+f)0 

0.20370E+00  -0.6]lllE+LiO  O.2037OE+0O 


STIFFENERS  IN  THE  x  DIRECTION  -  SX(I.J) 


-0.l06bfaE+00  0,l74b3r;*u0  -0 .  l0bf-6E+00 

0.106b6E+00  -0.l74b3f+00  U.l06f>6E+()0 


STIFFENERS  IN  THF  y  DIRECTION  -  SY(I»J) 


-0.1I636E+00  0»349U7r.  +  00  -0.  U6;*.6E+On 

0  •  11636E  +  00  -0.34907E+00  O.U636E+00 


FIGURE  A-IV-2  PROGRAM  PLTVIB:  TYPICAL  OUTPUT 


OlMfcNSION  Wl<3o>#*2(;i#3>  ,U<2#3>  *  V<2»3> 

Common  aq6)  #B( 71 
Integer  qs*o 
In=l 
1 0=2 

no  READ(IN#lub)  NDATa#IO()T 

ReAq ( In*  115)  MS#M,N<5,u»«S#Q 
RtAoUN#llO)  Al»A2»m#b2»H 
ReAu(IN#HO)  T1S*t2S»(,SS»ESiPRS 
HlAU(1N»H0)  T1K#12R»(,MR»lR*PRR 

105  FgRMAT(I4,I2) 

110  FoKMAT(5EH,4) 

115  FoRMAT(6I2J 
A(l)=  A1/A2 
A(2)=  H/Al 
A  t  3 ) =  H/A2 
A(4)=  fil/Al 
A(b)=  B1/A2 
A(6)=  t)l/b2 
A ( 7 ) -  H/Bl 
A(d)=  B2/A1 
A<9)  =  U2/A2 
A(1o)=H/B2 
A(11)=T2S/Tl5 
A(12)=T2R/T1R 
A(li)=TlR/TlS 
A ( I4 ) SGMR/GMS 
A(lh)=ER/LS 

A(lfa)=(1.0-PRS*PRS)/( \ ,U-pRk*PRR> 

CM=aMb*T 1S*A2*B2/ 1544 . 0 

CK=2.0293b6*ES*riS**3/(A2*B2*(1.0-PRS*PRS) ) 
CF=0.1591b493*SyRT(CK/CM) 

Ci=A(l)**3 

C^=A(6)**3 

Cs=2.0^A(l0)**3 

C4=A(9)**3 

Cb=A(ll)**3 

Cf>=Ad3>**3 

B(l)=(1.0+A(ll) )*(1.0+2.0*Ci+2.0*C2+4.0*Ci*C2) 
b(2)=C3*A(l4)*A( 13> 

U(3)=1.0+2.0*A(  i ) 
b(4)=C4*(1.0+2.0  +  A('>)  )  *A(  l2) 

B(5)=1.0+C5 

B(6)=2.0*C6*A(lb)*A(lfl) 

Do  210  K=(jS*0 
«HlTt(I0»2l5) 

NkITC( lOf 220 )  NUAJA 


FIGURE  A-IV-3  PROGRAM  BXVIB 


116 

117 


I. 

[  1?0 


2no 

■  201 

202 


IF-  ( I  OUT)  116.116*117 
WRITE l IQ »26Q) 

SnQ=SGN(K) 

Fq=FLOAT(K) 

Do  205  I=MS»H 
SnMsSGN(I) 

Fm=FL0AT(I> 

DO  205  J=NShNJ 
ShN=SGN(J) 

Fn=FLOAT(J) 

Call  stiff<sk#i#j.k> 

Call  massism# i»j#k) 
FrEg=CF*SgRT($K/Sm) 

SmM=CM*SM 

SkM=CK*SK 

IF(IOUT)  201»201»120 

tel  d » 1 ) =SNM*SNN*A ( 1 ) *a ( 6) 

Wid»2)  =  5jNM*A  ( 1 ) 
teid*3)  =  tel (1* 1) 
tel (2*1 )=  SNN*A(6) 
tel (2*2)=  1.00000 
tel (2*3)=  tel (2* 1 ) 

Wl (3*1 )=  Ml ( 1  *  1 ) 

W] (3*2)=  *1(1*2) 
tej  (3*3)=  Wl(l*D 
U  ( 1»2)  =“FM*A  (3)  /F'j 
U(lrl);  SNN*A ( 6 ) *U ( 1 • 2 ) 

U(  1. 3)=  UU.l) 

U(2*l)=  SNM*U(1»1) 

U(2r2)=  SNM*U(1*2) 

U(2,3)=  5NM*U(1»1) 

V<1,2)=-FN*A(10)/FQ 

V ( 1 # 1 ) =  SNM«A(1)«V(1*2) 

V ( 1 f  3) =  V(lfl) 

V(2fl)=  SNN*V ( 1*1) 

V(2f 2)=  SNN*V (1*2) 

Vx2*31=  LnN«V(1*1) 

Do  200  U=lf3 

DO  200  Jl=lf 3 

te'2  ( II » Jl  )=SNQ*W1  ( 1 1  #  Ji  ) 

Continue 
Go  TO  202 

WrITE( IOf 265)  If JfKfSKMfSMMfFREQ 
60  TO  205 

WrITE ( IOf 225)  I » Jf K,FRtQ 
Wk1TE( IOf 230)  SMMfSKM 
WRlTE(I0f235) 


FIGURE  A-IV-3  PROGRAM  BOXVIB  (Continued) 


WrITE ( I0#24Q>  ( (Ml <11, Jl) ,Jl=l»3) # I 1st. 3) 

WKlTE(I0»245) 

Wi<XTE(  10*240)  ( (Udlul)  »Jl=l*3>  »  1 1=1  *2) 
wp ITE ( 10*250) 

WHITE (10*240)  ( (V(Il,jl) »jl=l*3) *11=1,2) 

WRITE { 10*255) 

Wf,  lTE(  10*240)  ( (W2(I1,J1)  ,Jl=l»3)  ,  1 1=1  *3) 

20*5  continue 
210  Continue 

215  FoHMAT(iHl»7X»4bH|)YNARlC  ANALYSIS  OF  a  NIrE  CELL-'»Ox  STRUCTURE) 
2?o  Format (/»21x»9Hoata  CASt* i5) 

225  Format (/*3X,i3HM0uE  number  (*I2»ih»»i?,ih,*i2»ih) ,7x» 
HoHFREOUEnCy=»E12.5,3HHZ.  ) 

230  Format (/*3X*11Hm0dAL  nASS=*E12.5»3X*1rHM0dAL  STIFfNESS=»Ei2.5) 
235  Format  (//,22X*lbHi.;0r)AL  AMPLITUDES, /»  IrX,  lgHuPPER  COVER  SHEET  -* 
l7hWl<I,J) ) 

240  FoRMAT(10x»El2.b*2X,fl2.5,2x*El2.5) 

245  FoRMAT(/»2JX»13HRIBS  -  U(I»J)) 

250  FoRMT(/»23x*13HRlBS  -  V<i*J)) 

255  Fr  RMAT (/»l7x,27HLoWER  COVER  sheet  -  W?( I  * J) ) 

250  FcRMAT(/*bX,7H(M*N*'))  ,3X*i2hM0UAL  STIfF. *2X» 

HlHMOCAL  MASS»2X,14HHRE0uEuCY»  HZ»»/> 

255  FCRMAT(6X*lH(*U*lH»  *I1,1H*  *H»1H)  * 2X ,E12 . 5, 2X *El2 . 5* 2X »El2 . 5) 
Go  TO  100 
End 


subroutine  stiff<sk,i,j,K) 

Common  A(i6),B(7) 

Iu=i*i*d*J 

CiJ=FLOAT(IJ) 

Fp2=F  ( vj*  1  *  J) 

SKF=2.0*(A(l)*F(8,l,J)4A(h)*F(5»  I*D)-f^,0*A<l>*A(6)*F(4,I»  J)  )/F22 
SkF=1*u+SkF 

RkF=(F(3»I*K)42.0*A(1)*F(p»I*K)  )*AU0)4 
1(K10*  J»K)+2.0*A(6)4F(7*D»K)  )*A(3)*A(i2)**3 
SK=ClJ*(fa(5)*F22*SKF+f,(b)*RRF) 

return 

Ei.U 


FIGURE  A-IV-3  PROGRAM  BOXVIB  (Continued) 


Function  f(n»i»ui 
Common  a(i6)»b(7> 

Rl  =  (FLOAT(I)/PLOAT(J))*A(fl> 

K=Rin.0/Rl 

F=R*« 


return 

end 


SuRRO'JTlNt  MASS»(SM»  I » J» 
Common  aq6) #B(7) 

C 1 =FLOAT ( I > /FLOAT ( J ) 
Cp=FLOAT( J) /FLOAT (K) 

Ry=BU)+B(4)*Cl*Ci 

Sv,=B<1>+8<2>*RM*C2*C2 


K) 


RfcTURN 

Ll.D 


Function  sgnik> 

c=i.o 

Do  b  l=l'K 

c=-c 

Continue 

son=c 

rlTurn 

End 


FIGURE  A-IV-3  PROGRAM  BOXVIB  (Concluded) 
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DYNAMIC  ANALYSIS  OF  A  NTnE  CELL-BOX  STRUCTURE 
OATa  case  1 

MODE  NUMBER  C  i.  1.  1)  FREQUENCY:  0 .78729E+02HZ . 

modal  mass=  o.i4i26t-oi  modal  stiffness=  o . 34566E+04 


MODAL  AMPLITUDES 
UPPER  CO*Eo  SHEET  -W1(I#J> 
U.703b3£+U0  -U.113/5E+01 

-0.6l87bE+00  O.lt.OOOE  +  01  -0.61875E+00 

0  •  703S3E+-00  -0,lli7bL+0l  O.703a3E  +  0fl 


RIBS  -  UU»U> 

0,bl87bE-*-00  -0.1i)Ou0E+Ul 

-U.6lfl7bE400  O.lnOQOE+Ol 

RIBS  -  V<1*J> 

U.b6B7bE+00  -U.hilUUOt +U0 
-U.S687bE+00  O.bt.OuOL+UO 

LOWER  COVER  SHEfc.T  -  W2(I»J> 
-0.703b3E+00  U,1  i?,75fc>0l  -0 ,703^3E+n>> 

0.6lB7bE+00  -O.li.OUOC+Ol  0.61875E+on 
-u.703b3t*l>U  Q,lJ375L>iil  -0,703f3E+0 


0.61875E+00 

-0.6187bE+00 


0 ,56B7bE*0° 
-0.56875E+QR 


MOOt  kUMBER  (  1*  2*  1) 
MODAL  MASS=  0.1a697t-0l 


FREQUENCY:  0.14371E+03HZ. 
N.OUAL  STIFFNESS-*  0.l3bl4E+Os 


MOdAl  amplitudes 

UPPER  COvtR  SHEET 

-O.7O3u3E^00  -0.U375E+01  "2 ’ I? 0 

U.blb7bE+00  0.1U000E  +  01  2w5if7^t2n 
-0.703n3E+00  -0,ll37SC*0l  -0 ,703p3E+0f> 


RIBS  -  U(I*u» 
-U.61B75E+00  -O.lunOOE+Ol 
U.61B75E+0U  0.1()UOOE+Ol 


-0,6l875E+0^ 

0,6l875EtoO 


ribs  -  v<i*u> 

0.11375E+01  -O.lLliOOl -MJl  0,ll375E+0l 
0#I137bt+0I  -U.10000^01  0.11375E+01 


LOWER  C0VEh  SriEET  - 
0.703h3E+00  0.11375t+0l 

-U.blb7bE+00  -0.1flUU0E+0l 
U.70383E+00  0.ll37bt^0l 


W2(I*J) 

O^OSpSE+O*1 

-0.61875E+00 

0.703a3E>00 


FIGURE  A-IV-4  PROGRAM  BOXVIB:  TYPICAL  OUTPUT 
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c  This  program  evaluates  the  modal  mass  and  stiffness, 

£ _ FqEQUENCY  AND  MODE  shapes  OF  a  WEDGE  cTRUrTURE _ 

dimension  a<22>*F(6),*<2»3>»v<2Q»4> 

_ real  imuo) _ 

Integer  q»qs»qf»qmid2 

- In51 - - I - 

1 0=2 

_ iMt  l)S«2Hf>99 _ 

I K ( 2  >  = . 07888 

_ Im<3)  =  . 0395b _ 

IM (4 Is* 02408 

Im(5>  =  *01632 _ 

?oo  Read(  IN*  100)NDATA»MS»MF »ns»nf»qs»qf 

_ Rr  ad ( IN> lol)  Al»t31,B?,  alpha _ _ 

ALPHAS . 0087266*ALPHA 

Rf_A0(IN»101)T1S»Gm1S»E1S»PRiS _ 

read  < in# ioi ) T2S » gm2s  » E2S  *  RR2S 

ReAQ(IN»101)T1R>GmR>Er>PRR _ 

A( 1 ) =2 «* ALPHA* Al/02 

_ A(2)=B1/B2 _ 

A(3)SA1/B2 

_ A(4)=A1«B2 _ _ 

A(5)SA1/Bl 

A(6)=A(1)»»3 _ _ 

A(7jsElS*TiS**3/(l.-PRlS*PRlS)/12. 

_ A  ( B) =E2S*T2S»»3/ ( 1 *-PK2S*PR2S) /12* _ 

A  ( 9)  s£R*tlR**3/  ( 1  ,-PRr*PRR  >  /12 . 

_ A(1Q)sGM1S»T1S _ _ 

A(U)=AtlO)*A(4)/l544. 

A  (  12)sj,«+gM2S*T2S/A  do) _ 

A(13)=A(12)*(2.*A(2)+1.  )*Adl) 

_ A(14)=GMR»T1R/A<A0?«A(6)»A(11) _ 

A(l5J=24.3522T*A(7)/A(4) 

_ A(lb)  =  (l.+A(8>/A(7) )»A( 15) _ 

A( I7)=t ALPMA/1 *b7o795) **2 
A(l8)=A(9)/A(7)*A(l)*Aa7»*A(15) 

A(19)=2.*A(16)/A(2) 

A(20)s2*/A(17) 

- S7'2iJ  =  (r5.+a.'^RR) - 

_  A(22>=2.»A(20>  _ 

“DO  £03  MsMStMF 

SnMsSGN(M) _ _ 

FhsFLOAT(M) 

FM2sFM*FM 

- PlMoAs3.i4i59*PM - 

Dp  202  NsnS»nF  _ 

Sf>=S6N(N) 

Fn=FLOAT(N) 


FIGURE  A-IV-5  PROGRAM  WDGVIB 
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F,(2=Fr.*FN 

F  OKi<=FM/KU 

F  (  i )  =  (  (  FMoFi  j/A  (  d  )  +  A  ( 5 )  /FMoFf  j )  **2 )  *  A  <  l  c. ) 
F(2)  =  (  (FMoFfj/A(3>+A(3)/FMjFfJ**2)*A(  1^) 
F(3)=F<1)+F<2) 

L>('  20 1  Q=oS,  jF,^ 

0  1lj2=<0-i)/2 
Si  g»rSO  j(u) 

S'  vJv,i  =  ^Gfi( 0-11*32) 

F.  -FLOAT  ( .» ) 

F ,  2zF >,*F(i 
F(4)  =  (  (  A  ( ^0  ) 

F  (  S )  Ofa959*FM^  +  -  (  ?1 )  ( 22 )  *FQ2 

ho)=Mle)*(F(4)K  (^)  )/FU;.» 

Sr  MrF* '2 *F  2*  (F  (  3)  ♦Fft,)  ) 

S.  M=(  i.-o039t>4/F'  2)*  ,U4)*FN2/F5f2+A(  j3) 
F,.to=S  <K T  (  Sr-  1/S 
a  ( 1 »  1  )=Sfi..*A  (2) 

•i(l»2)  —  1* 

.»<1*3>=*<1*1> 

..(2*1)  =-w  (1*1) 

<i(2*2)  -“1  • 

<i  ( 2  *  3 )  -  a  (  *  *  1 ) 

0*  2u‘i  IK  =  1*20 

f-Fl0-»  r  ( ir* ) 

'i  z  0  «  0  S  *  R 

C...;4Z2  •  ♦  AlfH*\*S.*  4  ■  ♦*.(3) 

0  fVAz lu  ( Pi  -1  JA*  ( ) 

V  ( I.< » 1 )  -Si  ,.*COrt*S  -.f  A 

V(  Ik#2)=C  jM*Sv|F.<A 

V  (  1i<  »  3 )  =  V  ( I.<  *  1 ) 

<;  4  V  (  Irt » ■* )  =V  (  IH  *  2  ) 

».  !Tt(iO#2l5) 

«»,  I  TL  (  10*  220  ) ‘ISmT  . 

««>  1  Tl (  10 * 22s )•*•» > )»  ,*rF 
4i,  n ;. ( |o» «.3o ) s>' ••*  3)4 1 
...  1TL(  10*23o) 

•i.  iTuliOit^ulLOltJI*,:!*  i ) 

■»:  I  T »_  <  40*t43) 

ITL(  [0*. .40  )  ( (2.  J)  *v  zl*.') 

*,  ITU  10*2S0) 
u.  20‘.  1  =  1*20 
i\;Fi.O..  I  (  I ) 

^z  •  05*’ 

•>,  licit  I1.  >*t4i )  <*»  ( v  (’  *_  )  »J=1*4)  ) 

->  il  t(  *1  i  ■  i£ 

2'1?  Cv  *T1  .  E 

* o3  c.  < r ii i< •£. 


FIGURE  A-IV-5  PROGRAM  WDGVIB  (Continued) 
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\)  rv> 


inn 

101 

?lb 

??G 


GO  TU  200 
FORMAT ( 14,612) 

FORMAT ( 1H1 » 11X#  37rtDYNAMlC  ANALYSIS  OF  A  W^DGE  STRUCTURE) 

FORMAT  ( /3x  »  iIhMqOE  #  I2»  1H.  » 12, 1H.  ,  l2»  1H)  »5*»  IOHFRFWUE  CY  =  , 

.«3U1FcRMAT(/3X*11HMODaL  M,\sS=#Ei2.5»2X»16hM0DaL  STIFF  -jSS=f t.A2»5> 

r>3b  FoRi‘iAT(//20x»lbMMjU\L  AMPLITUDES# //16y »26H^PPER  CaVER  SHEET  -  m(1# 
1*1)#/) 

40  FV/r<MAT(10X»El2«5»2X,Ei2.5,2X#El2.5> 

41  FcRK,AT(F5.2»3X»4Ei2,5) 

^C;!^i7i:^HSruLLpUTLSH-  I,I.3!f;axaHSECT0„,6X.OHSECT0t<.h 

<;PL  ATE  lO#4X,artPLATE  Ni  # /1,3X » 1H1  #  i  IX »  1,2#  llX#  1H3#  1 )  X  #  lHu  » // > 

L.  U 


Function  sgnuo 

Czl.U 

00  b  I=1»K 

c=-c 

b  continue 
sgn=c 
return 
ti.D 


FIGURE  A-IV-5  PROGRAM  WDGVIB  (Concluded) 
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DYNAMIC  ANALYSIS  OF  A  fcEQfiE  STRUCTURE 


_  -UAiA-CaSE-  2flOQ _ 

MOHE  f  UMBER  (  2.  1,  1)  FREQUENCY^  0.2f.01lF+Q3>HZ. 

MODAL  MaSSs  &»_l9b22t*Q2  t  OQAl ..STIFF, NESS=  Q.S^lQEm 


MODAL  AMPLITUDES 

UPPER  CCtyER  SHEET  -  W(1»J>  ~ 

-O.IOOOOE+Ol  O.iuOOOE+Ul  -O.lOOnOE+ol 


_ LOWER  COVER  ShL£JL-=-JLL2±Ji _ 

_ O.lOQQOE+Ol  -0 .  lfjQUOE+Ql _ O.lQOnOE+Ql 


SECTOR  PLATES  “  V{I»J) 

_ SECTOR _ Sector _ sector _ sector 

R/A  PLATE  NO  PLATE  NO  PLATE  NO  PLaTE  NO 

_ 1 _ -2. _ 3 _ !L_ _ 


0.05  0,oia08E-02-0.6i808E-02  0.61808E-02-0.blQ08E-02 
0.1Q  0.23513E-Q1-0.23513E-Q1  0 .23513E-01-0.23513E-01 
0.15  0.48544E-01-0.4e544E-Gl  0,485R4E-01-0.4fi544E-0l 
0»20  _  0.76090E-01»0.76Qn0E-01  0.76090E“01-0.7ft090E-01 
0.25  b.lOOOlE+OO-O.loOnit+uO  0.10001E+00-0.1o001E+00 
0.^0  _  0.11413E+00-0.11413L+00  0 . 11413E+00-0 . H413E+00 
0.35  O.U327E+00-0.11327E+00  b,U327E*00-0.1i327E+00 
0.4{L  0»94052E"01~0«9uQs2E-01  0 . 94052E-01-0.94052E-01 
0.45  0.55627E-01-0.5bb27E-01  0.55627E-01-0.55627E-01 
0.50  0.537Q3E-0b»0.537n3F-06  Q.S3703E-06-0.5.37Q3E-06 
0.55  -0.67988E-01  b.b7988E-0l*0,679B8E-01  0.6798BE-01 
0 . 50  -0.14106E+00  0.1410BE+OOtQ.14108E+00  O.UlQBE+OO 
0.65  -0.21036E400  0.2l0 3bE*00v0.2i036E+00  0.21036E+00 
0.70  -Q.2bb31E+QQ  0.2hb3ir+00-0 .2bb31E+00  0.26b3lE+00 
0.75  -0.30002E+00  0.30002E+00.0.30002E+00  0 .30002E+00 
0.80  -0.3Q436E+00  0.3»43hH +00*0 . 3043bE+00  0^3n436E+00 
0.55  -0.27509E+00  0.27509E*00»0.27509E+00  0.27509E+00 
O.oo  -0.2n62£*ftQ_J,2Ub2EtM^211fe2E+OQ  Qt2llb2E+Q0 
0.06  -U . 11744E+00  0.11744E+00-0.11744E+00  0.1l744E+00 
LM _ -0.21481E-OS  0.2l4fllf-05<r0.2l461E-05  Q.2l4alE-Q5 


FIGURE  A-IV-6  PROGRAM  WDGVIB:  TYPICAL  OUTPUT 
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